=e alhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1978 


Stochastic separation of radar signals. 


Zach, Shlomo 


http://ndl.handle.net/10945/18553 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


(8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist Calhoun i named for Professor of Mathematics Guy K. Calhoun, NPS's first 


: \ KNOX appointed — and published — scholarly author. 
| inp 
LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 









% op ini 
eR ate 
9 (vein eekly baa Cu ° 
VES ie Ag rin’ vince RC Pan MAA Ng 
4 yi yh, is wy! re AS y af “Bhs ye . wi i Tye" Vie Oty WN nhs Vai ‘i iy 
' ie ne vw ,t ly rt . i eat 
t Hy \ : » i t i 4 e 


igi 
fh Mase ne Pity HNN, ut ie ‘yt 
a 

















* 
7 
ae vee NSO Nt 
ARUN ee 
: AS Ne ae a WAS tay H SS a te 
a CAVERN ay 4 ANS 


ARRAN uN ws 
hi 













CRS : 





















































































































































a 
SSS “ 
ts &e aye ang h: ae smu ee 
ry } ANGE A eit eta nt 1 “ vy wn 
: NR NS te Nia Ke Ma RONG Ae 
| ie hy AA NS NH oe EN CONRAD A Se 
' sta Ws A AN MIRA hs ae TANS Ne we wank ts I ih bet 
Mi a te Passes.) A Ay * Oy Bs AN nha ay Nae 
»Y; 5 ‘% * wh \ o% fe. 
, tt ar sy be hn 14 ity seein ; 10) 43 Des nas MN x eat sy chil die 
, ray es a nn Le 42 §, Mam ately a wv KS Nat aint tenants XS M4 NA ti ; Bahan seh 
4 ) ' © h?. 1 t ge 
” I ae ‘" ,Y yt t A iv ih } a me Ree: Ast aia Mt Bat aA Synch FS ve FA ie ie NONE fe 
¥ mre } « > hy Me " bud 1 ' A 
if ty 1 i A Y 1a is * WA Vi ‘ £8 it we iA as } H a a ey Ch cet A. 10. ft di +E: % SN Nah 
' CAYS Vaciy % tS Fo teh . Wf A vy w Whey Rs 4 ni ‘hy x ‘ a AIT Sy Hi, aN RCA Nt AGA ee * “Tg % : ch 
f i ? * a Ah pe a. 4 ba at f at j wk Ay The eve Ns ete tribe ‘ jASN Mas a sts Rate Ar be km. 
ty a ye 7 ' +> a i * re rh an i y ty 4 t NS Ha \ a aN isan SAN ae eS sie < eee a a. ee Tee 
‘ a ¥ ‘ 7, ah a | Ly 7 Cia *DeYy 2) 7 whrk 2 Psa! Ae SAMs, ©: 
4 Pia , are wee : 1 7 4 t Heters ye, "nh - ‘' int x é ‘>: hen age xy “i ote Nt = é oy SOA 
: ‘ NS a yet 4 ey wt 1, ths ph he ¥ sbinan ot ‘s i Ai "ey 4 i Be gn ne : SSRN Rah oe NN i ah nes are Shoal fs oe 
1 ty ae Lad ( Ina ri 8 Sth, " Yi; ss Ya “ 
} , Yess J ¥ Ne yg ie, we fot nt NePps Leith yee. ste aes et ee Ae %¥ ee wk 
Wo Pai, at a An Lt ey Od “ee i aon BARON $ ATHENS Sa aes < ues Sea Sh Mak x Abs RS ae 
4 I ° 2. re A : Fh A Mi er tA th yA hing Palen’ u X¢ < y yey ene EAS Shs sen otye a i £ ae ae a 
> a i A ee | RPL MR URS Ay WN tah RENE ST, Deak UREA ARAN GSI Write eT ries 
“4 ’ mee Ve ae Re 54% ni hvry & vi Wks tke Bikini ae AE ne ie SN yan A) ee WES Sa Maal Cae 
‘ ‘ a eee e 4 vs aihea! . 4 + St ene “" : yy a NS | CUE teh Mado Se Ne ah 
a ’ ‘ YW wy 4s ‘ 4 yh ’ a, " ws 1a et " * ‘ bath yt Siagten ss SiR, 
4 ’ ee Tra hh Sti qe AN A . ‘ x NOL eo Mee Gh OY Maa WEIN a Tey oy e Pt PTY Va Le EA RS 
i : > wry oi it ee Noe Vabaly ed W's hinge ist db py aD yt : rea vai " +3 Adc : ee Sa =< 
\ \ : Le ¢ \ me § ‘hh OGY a] ae Ny 5m 9 Vth bn yy Net mye habs xT 4 hee cn PSD Aa Tea Bae g 
‘ ‘. . \ \ nt 4 io “? i om Sy ¥ . J! ee ay ry AO NALA Ie he ‘ ah hs aN AREN itis tal Ke SS a age: nape Noh "alton, sh aon Cy en SAAS 
x a teu Le I MaRS TARY, 6 he i b eas Sh Thy rs yeNt ey RI ms ven PALE 
PR TOM NEN «08s SORE ae ER SERS SHR 
7 ¥ r 4 i 1 Ai me Way gw Wie Oyen ad Py ky NOSE SO, sesh y ee ane Le, ory, \ ah arty 
ot 4 a ; Ni ‘ LA ae roe nhs We FR Welty we arg AACS e+ Ae a tiekk inthis 
° % ‘ ‘ 4 wh a ‘ 4 p eR eke ake Y Ht ean ayy ‘yr bY AHR Pettit a 3a TAL Oe art ty 
- j Sg She i + i haan PAE N) ¢ <\ s ‘ . W ‘e ‘ - " sat he ™ Seay a Panay ees ayieen. ay SAN NSN bai pene 
’ ag i , : Vy ? au : Yan tol? 4 LOO eM Wi Phi SaaS) a het Bh ae . Ay B.A, Was a Thy ue WS). i 
a , { : ecu i b, Ape hh | St Ay ¥ “ ‘, : ‘AG x Pun yh eh ra. j ba Lh : ma aN ‘% Ee % ra % ie 
, t \ a ‘ol, oa ? 1A "ea wy % As 1 a KOCK) in 8 4 unk Ye \ che ih vie BBY ees Me tte Wt ee MATa 
Hr? 4 i I “a a i ‘s ; Mi reba 6 Sen AY is 1 hee: Fates as eh, v8, Rn WAR Ryle ee ‘ . hy cs, 
So ; ews tans WY Aaa A ako A wy Wana cis Pia, Fe rae A ah Rah Ricks ah PEE ant A ¥ ae at 
4 TA ae wr abe Ge SAAS Sky Pate ReMi oe hy 8 a ery he AME Bein PEW d bly a eA 1 Re al on 
' nr. : ef of “ent Bea PARA) ke te ANU a ye ON Aa Ravan shh s MS Sh cera Ray ieee aan Re 
: , i’ aa Wis "RRA Shea Maat yay, ” — PN Mise he Godan SoM Vrs ht ee, heme here Hs: 
atk Ay ae Ni: RY Say Wa eA RAS Ae ot i Cet -ey SWRASAR Eos to, Pose 
a_i \ ae: Be hey hee ca RNAS RSC SPAT: SN ane manarit Site SRERVA eon 
P h YY 6 a sie , Me > * * Meteo A ALD, Meee bp ie We Ay, wea the *, 
‘ : i Pat ayy PWM lek ull ay AeA 1 : Ft AMS Sgr SA. cant 
j aK Lar TS MV a my hey Py a Wat's Ny sahese Stn Sat erent {HART Sy, 4. oR eA tats, 
i Y 1 rhe weg ay War Ys ate hh Rea ae aN ae <(aert eat) Es Cea. Le a 
N, ; f ‘ “2 eth G a, te eee NY 4 a i ove fet ¢ \ He beh the tt Ay Ee aa) oe: 
\ ; : 4 «gy Ne (3 7 he ESA Os 8 se, x ry; a sehen. Ce eh 4B, MSC CaN ges BM re G eR Saatant ig dtl, ve, alk 
"Mle « ° "i ‘ - 5 Te ohh BS | fe WS Asbiteg tee FQ 
‘- % 1 4 wr oy) i eo UY * dal bk te pers anh Aa? ae pw Aad Aas ries ras eT Sty any Soh ro 
nk] , ee we . pIgt » “ ¢ < Ry. BO 945 St 5. Me ee an, "he c 
' * } ny be 5 vt aly Shh NEA : - ‘ Sete ; oe Ne Made gen : rhe Ae Ae iaon, WAw 
F ' ‘ » oh “) seh AY ee : 4 yy ; Fag me We ea yee A Wake ea Dhchs yak Ay pe 
<f § ' , acy Ls n i> i) er i 95 We ie b thd 
] Bh hay 8 ATS 5 4 SN ae a OAR shy Ce 4 
1 tt ; WETAN Se \ 
‘ ; ala dy Betws Sy St eS. RA 7 
, i, Ne et ee Ay ga ¥t ao sy 
: “i i 7) Mato.) Shh ae ey ra wd, 
vs , * ar ink EH, Se 
he. aha BS eh races OOH toy 
$ ; Kohshononch Woes hey 
. 4 5 ‘ me, FAGY FE 9 De. hae ay BIPRTRIN. 
| . SS SARS 
’ wie at Niet B+ eS RS is 
~ ae + % 4 £54 2A ae 
map ? 4 4 “yey 
{ F ‘ r) ote, % aT NE a,s% a a ae 
Ba ' 7 . 4 Us VA "RB. Nay ot Uy te 
feted, b] , + » 7" ry : lot d e aX qh A ee 4) hae) toe 
anal H x t > tate 4 iy a — res hy ‘* FM ~ i Bet ba 
to? ‘ + 4 » q: Lah Vs wt S >? e at 7 rn ey ale 
on 8 e : “A 4 ’, t + q Werle ae | tr teas wa 
J Ls "s . i 5 "1 ie r a's Nand 4 REE brah 
1 : f ‘ { . ri 4 ae ee xy 8a . ¥ Ree . = - a 
paw a: : ‘ j f y é ¥ % T 5 * a i & ; 5% Ett? » Saas 
a 4 > , ‘ ifs een’ : = i ' =m 4 LAL. 
aos ee i a he , > 7 0) 4, i ‘ A be ‘a P ate a Sa % a uA 
oa we - 1 m : d a is ad ame AY Seth iy q 
4 3 ‘ " a ‘ aa, =| * ae aie aU he Ae ods we f 3 
~ 3 ’ de i > ’ Mis ‘ ue is, " i ' “@ 4 tSghy ete nxt A. : 
< : be “x ; } ray . r “ ru Jy, 5 % 
‘b> 4 J rr %) ; , vA ow ; 4 Hays of *% a2 
. q ‘ Aa « £ y ly a ba) ¥ re * ; 
; ; ‘ i . % : 
' : : ' : Ly v eo ae al i et Se ate aks, 
a i. i, : a ‘ hee af gee ye e ™, Pe 
1 ¥ , ¥ 2% P ie aa ‘ a | 5 / 
ad ‘ 1 AA s os x ? OD) rs ‘ue ® "yy - 
7 7 , » b oe Pe a 4 4 f" 
ot F : 3 a i s & Z i ‘ ' + t% 
s a A _ 
op > Ww a - ey . 
. ’ ; ret . *, : a e.- | { ¢ ° 
# ; y ‘ : rr oP, - 
ew ‘ L hd . i oe f % worl ‘ fal, 3 ; 2 
a ’ , ' TAN »* ; ‘ ee ey mend 
a. , we r ail I P ‘ 4 $ ee aa & ty 
a of - ‘ A a! 4 ee bag 4 Ree” a Aa 
. 1 | aT | 
san ’ ? k vos . 4am a é Bia = bear eae Rw 
t 4 4 Wartime ? 1 +r) ~ ¢ Sere 
n ' 4 ta eer j J ¢ n 
- ‘ i W ‘ F . tad ite & BP Nya 
: : ty " a, ; ’ fet, 
~ ’ a « ey rd . ak Uy bets 4 
a ae : L Ss fe ot 
> ' = ’ os P , 4 = i HE alfesy pe ec J ae F t 
‘ * ” a ¢ <a wager Le f ‘ mars Caled 
F : . pan ‘ ; ‘4 to Dm | 1, it hy Oe ‘ eer Fos ay S mo eed . we bas = ag sg 
> Jqye Yas ‘+ ~ = i pen 
4 of dao, mn é ad on eyo ‘otigh, Bed as i wore lnye ting Vee J i See a pifmtnet gues 
ay . / Y > Pian ‘ r we aah a oe i oy Mies i % ey (re wary, eee ee Ming PS ges pe 
‘ , & mg a s Tots Bene AS ge | he ne tt ai ae Re mph a eee s tes vy gare 
r] 5 S : A, 2 ins y te Rig “ ies , ae chiar * oA Load ¥ eet Ee arre LADO es pla 
c. ' , 3 ' 9 ‘at 4 wy 7 a, .f ; ¢ suai UN se oe abst PRK ene ee ere “yy ae a nee, os 5 Sime em Me 
rs i! : ¥ y ; 4 F sae A he er Py hee “ug : oe, tt i tT v pd ay pha £4 £ shh ae =pr) 4 irs wes al ze 24, eae ri a 
9 é ‘4 ‘ £ ¥ x _ : a S ae bol sf Age = mS We Bx gs = ys i nm * hii ie ae bd 0 = 7 Hg 
ms H e ‘ j ‘, J 3 *, 4 ae ; A 5h tee e 5 ‘3 ' Cee din can x “Stee ven . eo * Hs 
F ’ ; 4 r ia’ 1% 1 te aw # “L be is Parse 
(Sa 1, Petes, 5 = orf 
* ; a } : ’ A " ce 5 a P| . ie . ue ip l= af m ” Gost a . ? se. ivr pr 4, peek (ne =v rae. iece teed sa: oS Ed? 
a . | Hoos ae Kopel endian”, ink ai hin ae dalle ie are borane th ao te 
' ee ee nal “Ch Ae ey . fy , es re) 
| i ' : " . aur", bre A (Q 5 1 f Pi "t 4 Ps 2 " A Be Was Se . se 1 7 ag Be salyapt mente) ce set sane ¥ ae at 
' » | t : a WE pb sf : oes - ta? H hp 6 
ae ' ‘ , 1 * * - * i," F 7 Bek Pp at AG i, . 4 een oy ¥ at area: i Ds oe we 
itp . : a 4 “4 et ‘al, pee we ‘ 7. Ce ee uy S2 Cpine matical 
+ : sos ° 4 ey PL tet, gl 4 : Asn’, 4) # ete0) ae ree — a . ay y 
‘ : i F 4 - re te * q , Wr gee ‘ = it “ i 4 se 
* P ‘ ‘ ? a ake o> - . . /* 5 ae, Ka Keine GR eH ye fia Ae ata Od wes Vir 
i AL ’ s i, ooh ' on ya ph gh At rn ted pal tn ie Ave 's Hr ve ree x 
4 1 ts i ; : Fi + ver * . 
4 f U + 4 Rr ; Oo ' + eae , t ae gt ni a ee BA re i + Mt « ley Lge ve Deas ton 
- f U = a” Me fon fy i ' Cr oe met Rie faery - ety gh 
~ aod a1 : a oF 1 j : Pg 4! 4 ‘id wi ee ogres yt ? fom 5 tian tive est aig 
. PE 1 gh & £ Kiba ' : ae Pyenik, ets: del a Behe ack 
2 2 Riz bs : ; Saw | 4 € ‘ y , : T Peasy ny rr ys alsa: $ an pea a (ee TT ate thats 
tan : ; t sth - P vg , 4 3 : eee TS ay ie Pri 45° * MF Pitan ue 
- i ; ¢ kh 4 ‘ 2 ' a A ee ee Gis pa deyetiyd f i ZO dy Pay? ! Ha rather cad RINE pease, 
= P 4 ; ; Fs a c 7) i 1 as arg ee , x o PY v. be WR Am yy ite tokens pa ¢ ood ik) w ¥ e oes Cee x or faine 
} ', ie ce H ee ee A i fora ping Polis Tar at: IAA pis A sede my ad a ee 00 ob ph i ita: 
a, ; : ‘a2 8 ' , : j Ar es of 4 in * HY ; te gf ny ' 4 Ku pes nite “on 8 gd etal ete te * ia Sig Reine Se SPE, 
| « a “ ; at oan ‘ 8 ‘ Maks j 
<> - : 5 7 4 a i u ' as i or t Au He fe oa “Oe 4 iota Peg oi tlie - pete” tsogtbprd Rsobenil Sapneng i re 
‘f ' ; am | te A : A - J 4 
a * ye a's f . Kes sa ° ¥ ' ” a pi) So 79 S, aie? "ty et el BD hat ff Fx) ‘ e. Figs r eat at E98 Ps dap 
2 , , PA ‘ beg t ; he ‘ yi ty: {+4 ay fom %, she a” AD iy abies ia tet pie C 
= ‘ " ‘# ae \ & :] ? Cee ' ’ 7 te ’ ? ; 4 rif % ‘ / iv. age hin Aine ry eee oe wr ath oe 
| oa 4 f 2 4 ‘ 3 4 ) hd ad “sn teh wh a FEAL od , CD | #4 : “ex 4 aiay & 8 diay i Paes 
ait x} oy ; ve vie é as \w faa dar eta if ag ot es fee foal ae ey EY 5 
ms . ' 1 r ‘ r , ft Pe me sy ' t ¢ fa Met . 4 rr) Poe leat F P ay 
‘ne is ' ' ‘ ~ - ft u i Ms! MRS ” APP oo - ad *# Pea ie rue 4 Rete ny aes ae ah Hines PP, Peg | oud “Peni % 
am » ‘ mo. : e : ; i yi if i Wet a : ¥ aa “i ” . ne : io ge rf 5, ah 3, ing 
# it ee - Nae ¥ . eae : ‘ : a) \ ess it, 4 ry . tT ae md, ‘Poke + ai wa*4 + At ip 8 Aah, re ay aa DF naan so Sts 
Ore yo 4 ‘- " mre ? ; Ls * fe ya of ee ‘ fed ue a, a - ¥ righ ‘if 4 Ene os a pe ItIF at aw Upaieagegen, rene & tear 
q : eae ‘ t ‘ ; ' tw ' 4 x“ OL Se fan pe 
4 Se, Wy Pat: } 6 ty * ae antl, if eae a pt it Meo F F “8 eye we ft ee Wy ; ge Be [Ps aA Way eet or 2 aaa ren 
- ‘ : J 4 ‘ of eT rea yee As 1 ie 4, wn ed ) at Pi tye f ‘ A Lit Peer 
4 j ' ; Ws t cy 4 1p ws “MUR seared wd lds ef Fra y ee va eee ms sbarrhg Sir Pepe ys, 
_ ‘ ’ ' ser’ t i rhe at We Ls AT 4 22a, * Sah 
J ry’ P rej bo a ‘ % any od ae es 7 Int ¢ Bit dy wren Sean) fi P fut 73 vo Mae sey, ayn 6, 45 Cb ay ts regh momar a 
heat Loe ‘ ‘ $ + 6b ta “ % a Og : ye ‘a ¢ a, Hey = 7 ‘ef PY 4 Fo Pacis #4 wir Pet poise,tog 
‘ ‘ P i 4 } 4 ' iE? e5ym, y re 0 ge he Ln a ope DLP al ; H sake, a. © hae ht dest 
: s : ; tee a ¢ t ' Wot ’ , Mair fr “4 wee Me ‘yee ee hy 4 nay» be i" Sam the, fel ghAL AI 
' a Y ' te Ns ie if Pie F - ba: hh pig ‘ be Mea nt mat Wed bia ih i ¥, f 98, sie gin. fh 
r 4 , . ts ro 03 Pliny, i Cer fos , et oi 
4 5 PU oT if ‘ ae Same 42 ; al Mae ah ry ge, Gpthad ok hares, ¥' 4 rks is ihe hd é' ; #440 16 39 EAS 
cy ; 4 ‘hu d ; Weusis.s , € iF Dei y as vitae. ap nite Va $4 Cia pe tet 5 do. fees: "2 L ve rt re A Wat aa Peng ‘Ow ly yg 
. J L 2 ye 1 ip 1” ] A i 
‘ : " : : ff ' ‘f sedan d's Lier “#4 i : er uy : *, vd ws YB on , Waits "f, Mes Fy /; ae, he # ve ! ete nai Ate a ns ! Pag fur 
ae bt fey ¥ en i ne er 7 erie Ay APT ae ne 
; ; : a rf ‘ oye f a rs $b he if < ae | at i jor i , 4 i 4 Pps wy ’ sdf 
er Pear. 4; gee ih a en Sey a ee ite oo gas a 2 hat ry rh SR AGTH Aye a Abs es aac ee A 
2 “e  » ? { Bh gi ae tae Pe Ms an ' F. Ady w rata ail Aer 3 ff ” i L cha ) Liha g rut Apa ete Ms 
” aia 1 7 pe a ; {i Ag eed 4h , Pag te designe tio a eo ¢F OR A se aBhA Je sof < tA". oes ae ney oy ee ote iP 
é ‘ { } PCs A Ese gy oe } pa Prk’ Piet 4 dhe Sf Tae ae oF 7 mies ePARS fis teaa gd Pet nh Se Ps Sarg lye, es ; 
; ae Su ye: P i Mba trysinte o Pt PE ie. Sa A MAN : ida. eae ag gai i Siawinag ea 
nity . 1 cé xt og a eee Saar Riek, ee fy bys slong sods te bef } nt ea ‘ 7 ; D3, pay vs, ge MS 6 las Santee 
: P cep ' te diy WP , t's Z A if, " ings taal © age ge Ob REO oe iy boa ; : oe : gan Pe eit Dy an actin aE te ay 
° : ' : aay ¢ pas ys ? no ee Bee Re SES! OMS hy And pe hog: me Mulviad 1) #1 é Cap pile a, Peston S aire Fi pe = 
F J, yO as rf  _<oae ; f, 5 2, * fy, “ Ck dios ahi fe if bteh y 1 nt Al dt Aa Dad vite ‘ Be ae When Pe Aryhsy i hist ie / Rt ety ie eee hed od Sees tie PSU cone p 
‘ es 4 ¥ F ’ { ’ ¢ | ae th doy uly Sy rity i? Ved Po ~ i rn rp ? o j : 
P Ni per dua > * has othye Cas : ett a tds pe ey ae 8? 5 8 tory ey aeeen - fae thew ts eer: ES doc 
, afc, ia uw tt Oe .. ental ‘ Jf baba it ah gaoy pest $b 4 he tt Ra CW é re us ee Phot t ie, fy a ra ee ret rane otigae > ate 
‘ ‘ : ‘ta RON 7 . * i ‘a fs Fa oAars y gts j oP iy od B lent, be Bt 7 ii Beit ether ty REE oo vs) Nepni Seis dake taf Seid ‘ ; pg: ‘ Pet : es 4 ; 
- = eet. a i. mr Te Fe! + tinh rite oy Oe if i? Fei fan i ny ea’ Ha be bt ae lay Ao Piet de re Age RE Peel aap ienaae : uP ids, potty 
, 4 a s "es ‘/ sp iy beag Hele jbiher vee if a any Vien fd ED " at a fb hay aa pws pf 
yt eat y 1 bebe it ’ hae ri ¥ : ? 
va! ee > nas fOp Pam oats, bp he Meg SPB Sahay ARR AES a or ele VA ee PEE Y, 
‘ : a ree eae ean aie i yepat ot Oe ae heed st) 4 fd ey yi, ny ime by sibel 43 yp ip Hi 2. %, had Lod yin: Ne ae oie 
f j $F y 2 i JS ! - SJ yy bbs ; . Laine el by ey Hy bal hee, Pr) Lane i Hy 
4 tf ch ,é af a, ' | Fee shea Fs 3. At ‘ pe ut Fs ia ihe iP a au 
' ; i) o) eey I Wn ‘ i DACA: 4 at nb aleds at | Rey py 
iy ee, | A De Fhe oe hor © Beg B78) dh kate fe i; wide i df b eS ae be é 
Aa! Med | att decd Olas Gt OPE OR Prien Rete 
b : a. Pay 4 ny ee are | aes * He, < $e baty? f we ev gt; an de ih, bey uA ae, 
Bi tia ‘ , akan Lt | Df HAP Site B20) a Nip: 
* 4 Ob: st act. Gi wy) Oy fica Me Je fey he PIVEN bys oxJb ile 
? J ‘ er day Wet ards Fs as yyy wg Be oe (54 fp ty 
eS a ; ihe 3" 1h An si! 1 age f saga fi 
- ; \ ra ae bist ye teat Nell Rr pr aes 
" ae) | ' . J A, vere % Age aes / ts a an A as y's islet 4 
" P fe 8} hi : ; 
i , : i 5 RD at bE 3 wm, rg i Lato ep Por ha wy TA : 
‘ ; Li . P* Ses : wae 7%, ‘s Sot Malh H ‘ots ae ‘ if bt ; ie wand bt 
a ‘ M7 eo .. ! a f ae ey "eles » ht ‘ i ; é - wan % pe 
* nf ; rs ft tr Gog Poe be LF rth 
" . 1 ,! a site ‘ i, ce AP wi yeh BAe Py! vd if pfs en ele ipintine ait 
i '} } ‘ “i? Ore PTS on dL vi kad i vy vented eed Haguey 
fae , : ‘ { f evar? rar ia ay fe SE 
, } . 5 ik t : -. ' te Ay ) vt atpey Ns bey ay 
5 a ia DOI Wr LE sc tyit 
: , : f ry re Mia Dar) sty ie 
Rca A aa a Is Fety Se 
. Sak adn i & yi Nyt 
‘ Tin am ey : i 
aa. , Si gh 1 Ue ie of aera aide 
ta 0 698 , Po reet at, Fee kh dtl 
‘ ¥ MM $i ty 4S Aa (? ‘' vr aly 
f ‘f yas ees Fad re, week ie 
' af ip 0g Cm ARS IO Se hy Al ef 
oo Lint tte Sef, fs nr, 
f , ‘ bien eet, EE at RAR ar gE 
oe | i ' : AY ¢ a} fl * sae $ | abe sf We aed fas nd £294 Ved 
5 ' gat pe ey BY ela 5 4 ie ke 
aie ‘ 4 ; PT AMPA ar) Pad NY Pa | Peay Buen! 
, 5 1 Us ve Layee ee a my Ce ad ea 
Lo tu a fase thy ‘ ' i ‘ ‘dt hey ; y Sgt he Ay S a 57 ihe ee if 
; . : fae i Fy Ad cht he Me Ne tia ght! rit ny 
A | ‘ha i eet Ve as Bs, ¢ Hog A. ? fay ay ba Tis a,’ ‘ ¢. Px "et vis 
’ ae oie , ' Dob 2 rf Ay; ee vis A ‘t, uf HC Cig : ny a 
4 Pe ‘ ee spies BE Hy OF ate hae) Pn Ory i 
1 i i Wy neue a fhe ries cn L Phen, hee AAS ABE) 
' a v's tet AL ' 4 i 8 , 
‘if ea Tiere UI ar i | ty fy t 
é : ooo) 4) 
ied . t, , 
ae ft, 


, r] Sey Lb 
ME had ee ahh t 
' a NU is hs and 
y *¢ 1, i] wh Wy t eh Baw Sa Md, 
elise & ti rei Fv: i aan dy } 
I 2 vs : + eri Teed | yee é ‘ t 4 
‘ ft ye i fo. - ey 
Pr anon A ‘ i AA ‘| ‘ 
RS 7 oe 
by me is ie F 
®, 


24% 6.37 Mel 
Be BC Aes Se ee ae See | 3 
of eg (fh 

' np Qari >= 

q is i . a4 s 

ae Tee? ae 

















ER 


C 


NAVAL POSTGRAQUATE SOHOUL 
t 


Monterey, Gallfornia 





STOCHASTIC SEPARATION OF RADAR SIGNALS 


by 


Shlomo (Welss) Zach 


December 1978 





Mies Ss Ad yiLSox: 


Approved for public release; distribution unlimited. 


1187412 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE 


"REPORT NUMBER 2. GOVT ACCESSION NO. 


4. TITLE (and Subtitle) 
stochastic Separation of Radar Signals 









READ INSTRUCTIONS» 
BEFORE COMPLETING FORM : 
ee 8 | 


S. TYPE OF REPORT & PERIOO COVERED 


Phep ww arhesis: 
Dec 1978 


6. PERFORMING ORG. REPORT NUMBER 















7. AUTHOR(s) 8. CONTRACT OR GRANT NUMBER(A) 


shiomo (Weiss) Zach 





10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMSERS 





9. PERFORMING ORGANIZATION NAME ANO AOORESS 










Naval Postgraduate School 
Monterey, California 93940 










12. REPORT OATE 


December 1978 


13. NUMBER OF PAGES 


224 


13. SECURITY CLASS. (of (hte répert) 
Unclassified 


ii. CONTROLLING OFFICE NAME ANO AOORESS 





Naval Postgraduate School 
Monterey, California 93940 






& MONITORING AGENCY NAME & ADORESS(!f different from Controlling Office) 










Se. CECLASSIFICATION/ DOWNGRADING ' 
SCHEDULE 






16. OISTRIGUTION STATEMENT (of thie Report) 





PeepEoved for public release; distribution unlimited. 


17. SCISTRIBUTION STATEMENT (of the astatract entered in Block 20, If different from Report) 


SUPPLEMENTARY NOTES 





19. KEY wOROS (Continua on raveree side if neceesarp and identify by dicck number) 


Beton, Radar, Chaff, Modeling unresolved targets 


20. ABSTRACT (Continue on reverse side if necessary and identify by biock manber) 


There has been a continuing problem of estimating a radar 
Signal when the noise and the signal have the same power spectra. 
iimtes iS particularly troublesome when one tries to resolve two 
close targets with a tracking radar. The purpose of this research 
as to show theoretically that there are practical ways to solve 
the problem. 


DO , oie 1473s cI TION OF | NOV 6815 OBSOLETE UNCLASSIFIED 


S/N 0 -O014- 660! a 
/ ae : SECURITY CLASSIFICATION OF THIS PAGE (Phen Data Entered) 





Ee TE a Ys ee 
ESCUMTPY CLASSIFICATION OF TwIS PAGE(When Nore Entered. 





m0. ABSTRACT Continued) 


Two approaches are introduced here. The first is to 
ea@@ea process filter to the radar, and the second is to 
design a new measuring technique with processing such that 
the signals will be separable. These two main approaches 
have led to several new or extended theories which were 
developed in the course of this work: 


emeae Probability filter. 

~ A modified Kalman filter (MKF). 

- A measurement modification technigue for monopulse radar. 

~ A new approach to the design of a monopulse tracking radar. 


Simulations were performed to check the two major 
@eee@eres (the M.K.F. and the probability filter). According 
to the simulations we can say that a tracking radar can be 
Mmeeertied in order to solve the problem of separating 
unresolved targets. 


DD Form. 1473 
ahs INCTASSTETE 
S/N 0102-914-6601 Z SECURITY CLASSIFICATION OF THIS PAGESMROR Date Entered) 





Approved for public release; distribution unlimited. 


Stochastic Separation of Radar Signals 
by 


SP LOnemuvelss)) Zach 
Lieutenant, Israeli Navy 
meee, Technion, Israel Institute of Technology, 1970 
lao taval P@seGraduate School, 1977 


Submitted in partial fulfillment of the 
requirements for the degree of 


DOCTOR OF PHILOSOPHY 
from the 


NAVAL POSTGRADUATE SCHOOL 
December 1978 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93940 


ABSTRACT 


There has been a continuing problem of estimating a 
radar signal when the noise and the signal have the same 
power spectra. This is particularly troublesome when one 
tries to resolve two cloSe targets with a tracking radar. 
The purpose of this research is to show theoretically that 
there are practical ways to solve the problem. 

Two approaches are introduced here. The first is to 
add a process filter to the radar, and the second is to 
design a new measuring technique with processing such that 
the Signals will be separable. These two main approaches 
have led to several new or extended theories which were 
developed in the course of this work: 

meeine Probability filter. 

- A modified Kalman filter (MKF). 

- A measurement modification technique for monopulse 

radar. 

- A new approach to the design of a monopulse tracking 

madar: 

Simulations were performed to check the two major 
Maeories (the M.K.F. and the probability filter). According 
to the simulations we can say that a tracking radar can be 
modified in order to solve the problem of Separating 


unresolved targets. 
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een CKODUCTLON 


A. SEPARATION OF SIGNALS 

In many cases we would like to separate two signals 
from a measurement which is a function of their combination, 
Such as their sum. One of the signals is usually "noise" 
and the other the "signal" or "information." Techniques of 
separating two signals if they have different power spectra 
are well known, for example, the Kalman filter or the Weiner 
filter. The case when the two signals have almost the same 
power spectra does not appear in the literature, mostly 
because each problem must be treated in a different way. 
The following are two basic approaches to this problem: 

(1) Use all of the statistical information about 
the signal to effect the separation. 
(11) Design measuring techniques such that the two 
Signals will be separable in some sense. 
The specific application to be addressed in this research is 
the tracking radar case where it 1s necessary to separate a 
target return from chaff. The sections which follow present 
a summary of the basic theory with reference to the litera- 
ture. Some new results have been achieved and are presented 
in the body of the text. They include: 
a) The separation of the two signals by means of a 
probability filter - Appendix D 


5b) A modified Kalmak filter approach - Appendix A. 
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c) A new measurement modification technique for mono- 
Pulse radar —- Section (IV.A). 
d) A new approach to a tracking radar design - Chapter IV. 
The following examples are given as part of this intro- 
duction to show that separation of two signals with the 
Same power spectra is plausible. 
The objective is to estimate the strength of a signal 
in the presence of noise when the noise and the signal have 
identical or almost identical power spectra. In those cases 
(where we can't effectively separate the Signals in the 
frequency domain) we will consider the processes of iden- 
tifying the signals in the probability or statistical 


domain. Consider the following: 


Baemple i-1 


Tt is desired to separate a Signal from noise when 

they both have a uniform probability density function 
(p.d.f.) of different widths. The amplitude of the noise 
and the signal are unknown and it is assumed that there is 
no correlation between them. The meaSurement is the sum 
of the signal and the noise. Thus, we assume that: 

1) The noise and the signal are statistically independent. 

Pee BOth have uniform distributions. 


3) The measured signal is given by: 


ae. 





It 1s desired to estimate the amplitude of the signal. 
Peomestatistics we know that Z has a p.d.f. which is given 
by the convolution between the two p.d.f., (see Fig. 1-1). 
By inspection of this figure we can see that the p.d.f. of 
Z cGontains information about the amplitude of the two sig- 
nals. Note that without prior information about Ay Ol A, we 
can only determine the two amplitudes, without being able 
to state which is signal and which is noise. The separa- 
tion technique would involve calculation of the p.d.f. of 


meee tsO, this filter can only give us the amplitude of 


the signal and not the whole signal. 


Example 1-2 
Consider the case where the signal is a Square wave 


between levels, (A, +A.) and As, where the switching 


time, T, 1S a random variable with given statistics. The 
noise is the same type of signal with levels (B, +B.) and 
B4- Assume that the amplitudes ratio is given: 
een 
(alee 3B) 
6) aes 


Figure 1-2 identifies the signal and the noise. Figure 1-3 
identifies the p.d.f. of the combined signal plus noise. 


With the p.d.f. for Z = S + N and the given ratios Ky and 


eZ 





AGS 


1/ A, 


1/ A, 








wees I-1. p-:d.f. of PVOmeTOndilc wand the p.d.f. of 
Eicmemstime (Refer to example 1-1). 
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CZ/A,4,B, 8B) CES MCM CCID mane Er. 
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«--B, —> in Bs - 











Pee 72, = er B. 
Po=P_(1-p,) at B.+ B, 
Z,=A,+Bo+A 
p3=(l-p \p a 2 2 f 
s/Py in NM se Bat B, 


Bye (Uap alalabin } 


Pico lae p.d.t. Of the signal plus noise. 


IE) 





K,, One can compute the amplitudes and even determine the 
Signal S uniquely as follows: 
Let us assume that A = Bl, then the estimated values 


of S denoted by S will be (refer to fig. 1-3): 


A alse 2 OE 2 


2 IL 2 
S = (1-2) 
A, +A, TL 3e Z. On ZA 


And for the case, By Stee: 


A alse TARO y 7, 


2 il 3 
S = r=) 
LY alae Z. or Zy 
ewes ¢ 
(1) In this example there is no assumption one way or 


the other about correlation between the signal 
and the noise. 

(11) The two cases (A, > By 6r A, < By) can be easily 
determined by checking the p.d.f., (the four levels) 
and comparing Ky and K.- 

(111) The levels can be computed easily by the four 


equations: 
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Je tee sae a = age oie) Get 1B 


2 2 iL le IL 
Z. = 
A oe A. ~ Bo fhe A, a Ba 
(i= 4) 
‘a 
A, + Bo + A, at A, > By 
Z. oo 
B. + A. > By ako A, 2 Ba 
2, = A. + AL + By + Bo 


And the cases A, > By Or A, < By can be determined 


by the two ratios: 


ee 


B,/B. = K 


(lv) The Signal can be determined even if the noise is 
the sum of a square wave plus random noise (for 
example, white Gaussian). In this case we have to 
introduce the statistics of the switching time 
T(t) and the random noise. 

In tne noise free case we can show the processor by a 


Simple block diagram given in Fig. 1-4. 
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In these two simple examples, it has been shown that 
in some cases Signals can be separated from the noise, or 
the parameters of the signal determined, even if they have 
the same power density spectrum. Thus, the idea of filtering 
a Signal in the probability domain seems plausible. 

A second possibility of separation to somehow change the 
measurement, seems reasonable. An example would be ina 
control system where we deSign our measurement such that 


the system will be observable (See next example). 


Example 1-3 
Consider the case where the measurement and the state 


equation 1S given by: 


phe) PAIS) es U (K) 


me 


II 


Z (RK) CX(K) + V(K) 

In the case where the system is not observable there is no 
way to estimate X unless we change something in the system. 
One way to do it is to change the matrix C, so that the 
system will be observable. Let's now consider a practical 
example related to the dissertation. We will consider the 
conical scan radar (see next chapter for details). Inoc.s. 
we often have the problem of cross talk between the two 
Channels (Az and FY). A Simplified block diagram of the 


radar is given in fig. 1-5. The system can be described by: 
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where all the vectors are two dimensional and the matrix C 


ae 


b al 


is given by: 


e(k) contains the errors in both Ay Slee tTeoe ana are 

the cross table factors in the measurement. 

Obviously where a = b = 1 the system is not observable and 

we can't eStimate the error (the target position). In the 

Case Of an observable system, but with a and b close to one, the 
regular Kalman filter estimator converges slowly. In order 

to improve the time constant of the estimator one suggestion 


mento Change Z to Z* by (see fig. 1-5 suggestion #1): 


Ze (K) = E - 2(K) 
where 
AL -a 
EB — 
eee He 


Zi 





This will lead to a new measured vector (Z'(K)) given 


by: 


Z'(K) = (1-ab) (e(K) + n(K)) 


This suggestion will solve the cross talk problem but 
1€ will increase the effects of the noise n. The problem 
will be solved only if a-b # 1. So another way must be 
found. This can be done quite easily by changing the modulator 
reference signals. The solution will be that C > I and the 
Statistics of the new measured noise will be almost the same 
as the given noise and this is indeed an improvment of the 
System performances (see fig. 1-5 suggestion #2). The 
details are beyond the scope of Chiesa ciate. 

The situation of having the signal and noise with 
approximately the same power spectra occurs in radar. For 
example, we may have: 

(1) Target plus clutter. 
(li) Target plus target in the same resolution cell. 

The work which follows will concentrate on the tracking 
radar case. We will treat the problem with two major 
approaches: 

(1) Assume that the radar receiver is given and we have 
to process the output signal. This will give us the 
possibility to try to work the problem using the 


first proposed technique. 


Ze 





(11) Assume that we can change the receiver and also 
Make some changes in the antenna such that we will 
get a new measured Signal from the receiver which 
Will allow us to separate the two signals coming 
from two different sources (targets). This is the 
Second proposed technigue. 
Because the second approach will give us the most power- 
ful results, we will concentrate mainly on this approach, 
but we will include all of the information necessary for the 


first approach. 


B. APPLICATION TO TRACKING RADAR 
In radar tracking we can have the real target and the 
false target (chaff, clutter) which are in the same resolution 


cell. (AR<5 At, Aa<O. ap: where At = pulse length, § = beam- 


3db 
width of the antenna.) Because of the false target, we 
cannot measure the position of the real target 

(angles and range). For tracking, we must be able to 
measure the error in range, azimuth and elevation where the 
error equals the deviation of the actual target position 
from the estimated position as given by the midpoint of the 
Mesolution cell. In order to accomplish this, the radar has 
two (at least) windows for each coordinate. The "weight" of 
the target in each window is measured by the difference 
between the normalized weights of each window for each 
coordinate, which is used in calculation to determine 


target position. 
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Poample 1-4 


In range we have two gates, the early and late gates. 
By measuring the weight of each gate we compute the error 


in range as: 


ee = “ae awe 
ay a + V 
EG LG 
l=) 
Ro = ¢ NCI 2 
where: 
(1) AR - the computed error in range. 
oe) «6«2At —- the width of the range gate. 
a) the position of the range gate. 
oy) C Z—veloarty OGL light. 
e 
O 
aw) x5 = ie Vette ede 
eo daet 
O 
Etat 
ev. = f Wiel che 
te 
O 
(Vil) V(t) - the received signal. 


The time representation of the range gates, received 


Signal (video}, Vig and Vig! Pome Veneelin ho 1-5. 
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The same type of equation holds for the angle. In 
this case we have two or more antennas for azimuth, and 
two antennas for elevation. In this technique we estimate 
the weight of the real target in each window, and compute 
the error to close the tracking loop and try to reduce the 
errors to zero. When we have two targets in the same reso- 
lution cell, the computed error would depend upon the 
"center of gravity" (radar center) of the radar return of 
the two targets. Because the two targets have the same 
(or almost the same) power density spectra, we cannot dis- 
tinguish between them by spectral analysis; thus we have 
the problem discussed in Section I.A. 

The object of this research 1s to eStimate the true 
position of the target when the false target is near the 


true target. In the research -we avoly the following 


Mmeactical constraints: 


(1) An eStimator which can be used cn small missiles. 
We need simple algorithms so that the memory 
Size and the amount of computation will be 
reasonable. 

(11) The size of the antenna is fixed, and we are 
restricted to an antenna of reasonable structure. 
For example we can use a monopulse antenna with 
changes (referred to regular monopulse) but we 


cannot use an antenna with more than four outputs. 


AG 





fei 1 ) 


(iv) 


We will not search for a solution in the r.f. 
range because for each r.f. it is possible to 
have counter measure chaff adapted to that fre- 
quency or frequencies. 

We will restrict ourselves to surface targets 
which imply no reasonable Doppler shift between 
the two targets. In the problem of air targets 
the difference in velocity between the target and 
the chaff is distinguishable and they can be 
separated easily by a moving target indicator 
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Ii. BACKGROUND THEORY 


A. INTRODUCTION TO TRACKING RADAR 

Before we establish the problem and our solution 
we present a short summary of pertinent tracking radar 
theory. 

The function of tracking radar is to select a particular 
target and follow its course in range, angle and sometimes 


frequency (a direct measurement of velocity) coordinates. 


B. TRACKING BY RADAR 

Tracking systems can be achieved by two different 
Peenwhniques. 

Track-while-scan (TWS) 

This method of tracking basically uses data from search 
radar. The idea is to take a sample measurement of the 
target each time that the antenna is pointed to the tracked 
Meeget. This tracking class is not of interest in this 
Seady. 

Continuous Tracking Radar 

In this method of tracking, the antenna is always 
directed at the target due to function of a control system. 
There are several methods of tracking (i.e., method of 
generating the error signals to close the antenna position 
and the range loops). The important methods are: 

(1) Gonical scan (C.S.) 
1 ) Lobe switching (L.S.) 


(ii) Monopulse. 
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Those methods are distinguished by their difference in 
angle tracking. The principle of range tracking is almost 
the same for all the radars. Since the most accurate and 
interesting case for us is the monopulse, we will concen- 
trate on it. However the theory of the other techniques 
is summarized. 
meeangle Tracking 

peeee@onical Sean (C.S.). In C.S. syStems, angles 
are measured by a single antenna, whose radiation pattern 
rotates periodically about a certain axis. When the tar- 
get is in the axis direction, the radar will have a constant 
return signal. The coordinates of the off axis target is 
determined on the basis of comparison to the envelope of 
received signals with a reference Signal. For this reason 
we can not measure the error in one pulse but we have to 
Walt at least one period of antenna scan to get the informa- 
tion. A simplified block diagram of the C.S. system is 
given in Fig. 2-2.1. The antenna A, scans in space with 
angular frequency, %. Two references for the phase detec- 
tor are generated by the antenna (AZ. ref. and EL. ref.). 
The position of the antenna is controlled by the control 
system. Assume a target at position T. Because the target 
is assumed to be off the axis of rotation, the echo signal 
will be modulated at frequency 2 (the C.S. frequency). The 
amplitude of the signal depends on the position of the 


target with respect to axis of rotation, and the phase will 
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depend on the direction of the angle between the target 

and the rotation axis. When the antenna is on the target 
we will get zero modulation at the receiving signal (see 
Fig. 2-2.2). The modulation pulses pass through a receiver 
Wee tS controlled by an AGC (automatic gain control). 

We will see later on why the AGC is needed. The signal 
then passes through a "box-car" circuit which changes the 
amplitude modulated pulses into a signal close to sinusoid. 
This signal passes through two-phase detectors (for Az. and 
EL.). The outputs from the phase detectors are the errors 
in EL and Az, and by using these errors we can close the 
loop on a control system which will move the antenna to 
reduce the error to zero, sco that the antenna "looks" towards 
the target. We can write the amplitude of the signal at 


the output of the receiver as: 


St (t,9,) = S,°Gav(6,,) [l+m(o,.) cos (2t+o ) | (2.2-1) 
where 
S'(t, 9,7) = envelope of the amplitudes of a 
Sequence of radar return pulses. 
So = the amplitude of the signal (dependent 


on the RCS (Radar Cross Section) of 
the target, the range, the gain of 


the receiver, etc.) 
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Gav (o,) = the average gain of the antenna, in 
the direction of the target, which 
can be determined approximately by 


G(8 0 - On) 4 G(68>+ 8) 


Gav (8,,) = 5 


(22-2) 
to a first order of approximation, 
this average gain is constant, l.e., 


independent of Ome 


Oo = The position angle of the target. 

Q = Conical scan frequency of the antenna. 

Pe = Phase of signal relative to C.S. 
frequency. 

m(6,,) = The modulation index which 1S approxi- 


mately equal to: 


S (Come) = So ee | i) 
m(6,,) nex (Do 3) 
Eh 


G(8>- on) = G(8. +o.) 

Gio. -6.) + G(6o +06_) 

G(8. Sn r G(60 + ¢,) 

where (see Fig. 2.2-2) 

S - is the maximum received signal 


max 


= So (GS 0-om) 1: G(60+9,,)) - 


oe: 





me - 1s the minimum received signal 


aiacral 
= Sy (G(8 0-o mn) - G(8 +8) ) 
S +S, 
_ max uquibsa: 
S Sy Se ee ed 
av 2 


Because we would like to have the same amount of error for 
the same error angle, independent of the target and range, 
we have to measure the modulation index. The method which 
is usually used to achieve this objective is the AGC (autc- 
fees Gain control). The purpose of the AGC is to f1x ae 
omen that er will be a constant. This can be accomplished 


by controlling the gain of the receiver. The signal can be 


represented by: 
' a2 ) —4 
S (t,9,,) CL sp m(8,)cos(LQt+o.) J (2.2-4) 
where C is a known constant, independent of the target RCS 
and other parameters. 
If we pass the signal through a highpass 
filter, we will get (after the Box-Car): 
S(t) 9n) = Cm (9,,) cos (Lt+o _) (22) 
Now we have a Signal with its amplitude pro- 


Pecelonal to mS), and phase equal to oe With both of these 


we can determine the error signal after the phase detector: 
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AAZ = Kem(G,,) coso. 
(22—6)) 


AEL = K-m(6,,) sing. 


Because m(9,) depends on the antenna parameters, 
the error measured by the system will not be a linear func- 
mien of Opn The error signals can be determined linearly 
only for small displacement angles. 

tose eovercening (L.S.) The method of L.S. is 
achieved by switching the antenna beam between four posi- 
tions (two positions per direction). We will describe the 
method in one dimension since the two dimensions case is a 
Straightforward extension of the one-dimensional situation. 
Refer to Fig. 2.2-3. Because we have two antennas and we 
Switch between them, the receiver output belongs to antenna 


af (S,) and to antenna 2 (S As in the C.S. system because 


A 


of AGC, Si 5 2 S. equal a constant and the error can be 
determined by the difference between Sy and Ss Ae, 
error (@,,) K(S, (8,) = So (Sq) | (2.2-7) 


If we have four antennas, by AGC, we keep 


Sy eee 5 + 5S, equal to a constant, so then: 


2 3 4 
A = K[(S4 ~ S5) - (S, + S,)1 
(2.2-8) 
So K{(S) ~ S 3) - (S. ste S4)] 


Sis: 
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c. Monopulse Method. In the two previous methods 
(C.S. and L.S.), the measurement of angular error in two 
coordinates (AZ and EL) requires that a minimum of three 
pulses be processed. In practice, the minimum number of 
pulses in L.S. usually requires four pulses (one for each 
antenna). C.S. usually requires much more than four (20 
pulses is a "good" number - which is dependent on the fre- 
guency of the C.S. and the P.R.F.). The requirement for 
both of them is that during the scanning or switching the 
Signal amplitude does not contain additional modulation 
components other than the modulation given by the antenna. 
If the amplitude contains additional modulation components, 
the tracking accuracy might be degraded, especially if the 
frequency components of the fluctuations are close to the 
C.S. frequency or the L.S. rate. These functions* can be 
caused, for example, by a fluctuating target cross-section. 
In order to avoid this problem, we prefer to base our 
measurements on one pulSe rather than many. When a Signal 
1S received by two different antennas, the difference between 
the two Signals might be one of phase and/or amplitude. 
This is the premise of the monopulse radar, i.e., use the 
difference between signals arriving at the same time by the 
antenna. The name monopulse is used to describe those tech- 
niques which derive angle-error information on the basis of 
a Single pulse. (Note: In some cases monopulse is called 


"Simultaneous lobing.") 


a7 





Two major techniues are used in monopulse: 
ea.) Measuring the relative phase of the carrier of 
the I.F. between the pulses coming from each 
aleetina - 
(11) Measuring the relative amplitude of the pulse 
received in each beam. 
Fig. 2.2-4 presents a block diagram of general monopulse 
(1.e., amplitude - phase type), for one coordinate. We 
will see later that by zeroing the distances between the 
antennas, i.e., d= 0, the system reduces to amplitude 
comparison monopulse, and by zeroing the tilted angle, 
1.e., om = 0, it reduces to phase-comparison monopulse. 
The extension into two dimenSions is straightforward. In 
most of the cases the sum and the difference are taken waeh 


the r.f. Signal (right. after the antenna), after which 


mmenwe are two receivers (one for sum channel and one for 


the difference channel). This is because of practical 
reasons: 
(1) Only the r.f. channel receiver must be 
maentical . 


(11) For two-dimensions we Save one receiver (3 versus 4). 
In the analysis which follows we will show the 
Sum and difference taken in the 1.f. This change does not 
make any difference for monopulse radar but it 1s necessary 
for the complex angle (CA) method. The signals received 


from the target at antennas AL and A, are fed to the receiver 
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and then the sum ()) and the difference (A) are forming. 


Thus we can write the following cos term of yw 


Meret) 


A(S,7¢) 


where: 


IF* 


$1 +8, = KyS, (6, (8,) cos (wi pt + ) 
-~ G,(5,,)c0S (wept - >) ] (2.2-9) 
S,-S, = KS, (G, (8p) cos (wrpt + o> ¥) 


- Gy (On) COS (Wot -o- v4? ] 


= Target's position. 
= Fixed antenna offset angle. 


= Gains of sum and difference channels 


respectively. 


= The average constant signal in the sum 


Channel determined by the AGC. 


= The phase shift between the sum channel 
and the difference channel caused by 


Ene ClECULES < 


= The phase shift in antenna Ay and A, 
caused by the displacement of target by 


angle Onn relative to the bore sight 
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@epeectien (see Fig. 2:2-5), and is 


given by: 


1? S5 = el ais Seng ye 1c) b/elaaligaer| 


Single radar return. 


1 R +d sin oan 


aw 
I 


R 


9 R- d sin On (222=103 


4td 
A 





ell 7 
o = Swaeasl R,) = Sin Om 


2d = Distance between phase centers of 


antenna AL and A, (see Fig. 2.2-5). 


G = Radiation patterns of antenna Ay and 
A, respectively; for most of the cases 
they have the same shape only shifted, 


SO we can write: 


| 
Q) 
oo 
CD 

i 
CD 
ee” 


G, (8) 
C22 lee) 


G. (8) 


it 

Q 
@ 
+ 
@D 


where G is a Standardized radiation pattern, and §$ iS an 
Sieeittary direction. 
Let us assume that the gains in the sum and 


the difference are the same (not a necesSary condition): 


= ie = kK (2.2-12) 
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Sco Se . 
Figure 2.2-5. 'tlavefront Pela Gnomon i ps 
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By means of the AGC we fix the average sum channel to 
memein constant, as in C.S., we can assume that after the 


AGC we have: 


Ko So = constant (2213) 
Thus, for any level of signal we will have the same output 
error characteristics. After the phase detector we have 
an output which is the average of the multiple signal 
(} - A). It is shown later on that this quantity is pro- 
portional to the error between the direction of the target 


mative to the antenna. 


error (0,) = K' - Average {}(8,,t)-A(@,,,t) } 
= a [G7(,-8n) - G7(6.+6,,) ]cosy,, 22s 
+ [20G(8 -6,) °G(6,+6,,) sin2o]siny. 
where: 


K' is a constant determined by the AGC. 
a' is a known constant. 


From the last equation we can see that the 
meeor 15 only a function of on and equals the deviation 
from the bore sight (the rest are known constants determined 


Only by the system and not by the target). So from 
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ae 2.2-14 one can determine a Tomcsiioldry, Ghe computations, 
we can choose the parameters ae d and Te In the two major 
types of monopulses, we choose these parameters: 
(1) Amplitude comparison monopulse: 
Choose Uy = 0O0orde= 0. Practically, we choose both 
of them equal to zero. We would get, for the case 


we ="O- and d = 0: 


O 
error (9,) = [G7(8_-8_) - G*(8 +6,)] OPiS) 
AL Onl L Ove. t : 
(11) Phase comparison monopulse: 
Choose ve = jee, Or oe = 0*,. Practically, we choose 


both. Then, for the case oe = 0 and vs = 1/2: 


2 


error (o.,) = 2aG (6,,) sin UB 
= 206° (6,,) sin(ysin@,,) eo) 
.; Q 8STd/A 


From the last equation we see that oo can be 


determined. Near the origin (9,>*> 0). Then we 


A 


can assume: 





* 
Assuming G(8) = G(-8@). 
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aoe = 
G(8,) eG Qa i) 


and 
°Ss] oe 9 
Sin(y -Sin on) Y 7 
so that: 
On = B -error (8,,) (222153) 
where 8 = known constant, where the error is the 


output after the phase detector 


_ 2) 
(se-— 1/aG. ae 


When we wish to determine the position in two 
coordinates, theoretically we need only the addition of one 
antenna; practically, four antennas are used. The theory 
is almost the same; the only difference is that now we 


have three channels: 


( ) Sum channel = > 

(3) Dimtenence channel in AZ = Anz" 

meee Ditference channel in EL = EL’ 

where (refer to Fig. 2.2-6): 

we see Coes. eS 
: ere 2 ee oe Ohad od 
A = = = 
AZ (Sy oF S3) (S. + S,) (2.2-19) 
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oem ecu Tye locabion of the antennas 
in two-dimensional monopulse. 
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See serene CUtpUut Erom antenna 1. 


Bee Range Tracking 
The most widely used technique for tracking in range 
is based on two range gates. ‘There are no special tech- 
niques as for the method of angle tracking. Information about 
the range is taken in most of the cases from the sum 
channel in the monopulse radar and from the incoming 
Prmeses in C.S. or L.S. We will briefly cover the prin- 
mere Of automatic range tracking. A block diagram for 
Baeerange loop tracker is given by Fig. 2.2-7 and the 
time representations of the signals are given in Fig. 2.2-8. 
We have two major signals in the loop: 
(1) The video pulses - S(t) 
They contain information about the range position 
of the target. 
(11) The range estimate - R. 
The main idea is to open two gates, Gi’ So before and 
after the estimated time of the center of arrival of the 
pulse (t = <, R = target range, C = velocity of the light). 
One is the early gate (gi), and the other is the last=- gate 


(g,). The portion of the signal contained in the early gate 


is subtracted from the portion of the signal contained in 
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the late gate after integration. The error is calculated 
after the end of the late gate. The magnitude of the error 
signal (4R) is a measure of the difference between the cen- 
ter of gravity of the video signal, S(t), and the estimated 
range, which is the center of the two gates. The error is 
fed into a control system to estimate the range. The output 
of the control loop is the estimated range, R. To convert 
the range into time (to generate the gate's pulse) we uSe a 
converter which has the time of the transmitting pulse as 
a refernce input. The outputs of the converter are: 
(1) The two range gates (gy and G5) 
(11) A reset pulse to reset the integrator before the 
mange Gates. 
Because we want to have the same output error, AR 
for all the targets, independent of the target, we must 
normalize the pulses, and this is done by AGC, which has 


already been applied in angle tracking. 


one 





Pole ENT OF STHE PROBLEM 


fee RODUCTION. INTERFERENCE CAUSED BY TWO TARGETS 
Let's consider the problem of tracking a target when 


a false target is in the same resolution cell, il.e., 


0 
R 3db 


AR < <At-e/2 


where: 


Gi = the beamwidth of the antenna (in radians) 
ZAt = pulse length 

R = range of the target 

AL = the difference between the true target 


and the false target perpendicular to 


the line of sight of the antenna 


AR = the difference between the true target 
and the false target in range (parallel 


to the line of sight) 


o. =eehiesve locity Of ache. 


For Simplicity, let's consider the two-dimensional 
case only (range and elevation). The principle in three 


dimensions is the same but more complicated, and the results 


ae 





can be obtained by a straightforward extrapolation of the 
two-dimensional case. We will consider the case of point 
targets for developing the theory, but we will implement 


the results when the targets have complex structure. 


EB. INTERFERENCE IN ANGLE CREATED FROM TWO TARGETS 
Consider the case of two point targets (or spheres) 


in the same resolution cell of the radar with spacing 


AL Cie 
Pie (l.e., a < 9 3 Gp! a\iee << 5 


As the relative path lengths between the radar antenna 





Mmeetotee. Velocity Of Tight). 


and the two sources vary (l.e., Ry and Ro), the two sig- 
nals will alternately add and subtract, and so the amplitude 
and the phase of arrived Signal will change. Although such 
a Simple Situation (two points targets) may be fictitious, 
it will illustrate the main behavior. The relative ampli- 
tude between the RCS of the two targets iS assumed to be a 


constant "a" and the relative phase difference also constant, 
a. The difference in phase is due to difference in range 
Mmeor to reflecting properties. The relative angular 

error (AG/6,) 1S given by [23]: 


2 
Ag aio a cos @ Bao 


=< 
Cel = 


2D lL + a? + 2a cos a 


The position of the stronger target corresponds to 
AS = 0, while the smaller target position is at A6/6, = l. 


(See Fig. 3.2-1.) The position of the tracking system depends 


2) 





can be obtained by a straightforward extrapolation of the 
two-dimensional case. We will consider the case of point 
targets for developing the theory, but we will implement 


the results when the targets have complex Structure. 


tel PYTERPERENCE IN ANGLE CREATED FROM TWO TARGETS 
Consider the case of two point targets (or spheres) 
in the same resolution cell of the radar with spacing 


AL erie 
pe (1.€C., aoa 83 ab! AR < xi C 1s BNe Velocity Of light) . 





As the relative path lengths between the radar antenna 
and the two sources vary (1.e., Ry eMichee gene (Ewes S1c— 
nals will alternately add and subtract, and so the amplitude 
and the phase of arrived Signal will change. Although such 
a Simple situation (two points targets) may be fictitious, 
it will illustrate the main behavior. The relative ampli- 
tude between the RCS of the two targets is assumed to be a 
constant "a" and the relative phase difference also constant, 
a. The difference in phase is due to difference in range 
(AR) or to reflecting properties. The relative angular 
eeror (A6/6,) is given by [23]: 
Ag a’ ia eco s 1c 
= (3.2-1) 


OD 1 + Be + 2a cos a 


The position of the stronger target corresponds to 
A$ = 0, while the smaller target position is at A6/S, = 1. 


meee Fig. 3.2-1.) The position of the tracking system depends 


a2 





on the relative phase a and the ratio a. One can show 
Biae ©Or 0 < a < 1 and 0 < a < 2m - the values of A6/A8, 


will be between 


= 00 < — a 0.5 (3 2— 1) 


When the echo Signals are in phase (a = Q), 
the error reduces to cant which corresponds to the so- 
called "center of gravity" of the two targets. 

Now, when we have a complex target, i.e., the 
ratio a is a random variable which changes from pulse to 
pulse. Aliso a is random variable which changes from pulse 
to pulse Taking these into the given servo loop of the 
radar system, one can compute the statistics of the random 
variable A@, and this can give us a good approximation of 
the error (for example, mean and variance) caused by the 
addition of two targets in the same resolution cell. How- 
ever, the simple case results are sufficient for our problen, 
i1.e., to show that the result of two target returns in the 
Same resolution cell causes an error in the estimation of 
the position of the true target which depends on the signal 
ratio between the two targets. When the two targets have 
a complex structure rather than a point, the results are 
Much more complicated and we have to take into account the 


Circuits involved and the statistics of the target's returns. 





C. INTERFERENCE IN RANGE CAUSED BY TWO TARGETS 

In range measurements we have a similar effect, called 
glint, i.e., we will get an error range depending on the 
difference of range between the targets and the phase 


difference. We can write the measured error in range as: 


ae 
- ie AT 


AR = K[ jesce) des S (lei che! (3.3-1) 
COT At te 


where 


S(t) = S,(t) + S,(t) 
NT = the width of the range gate 
C5 = the estimated range of the target (l.e., 


the position of the range gate). 


The behavior in range is strongly influenced by the 
radar parameters and the circuits used for measurements 
(for example, AT), but a "good" estimation of the behavior 
can be taken as when the two targets are totally in the 
range gate ({l1.e., the two targets are very close in range, 
Or there is a large range gate compared to At), so that the 
range gate will be at the central gravity of the two targets. 
Also we will assume that the video pulses are rectangular. 
Merer tO Fig. 3.2-1; the center of gravity of the two 


pulses will be: 


i 


2 Se) ce 7 | s(t) dt 


ay 


Teg 
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S(t) 





Amplitudes of the two signals, separately. 


ce = 
Sy = The two signals combined in phase. 

S. = fhe two signals for combined out of phase. 

S23 = Received signals for the general case. 

— PReseencer Gl Gravity of the target. 

ie 5 The center of gravity of the two combined pulses. 


Figure 3.2-2. Range errors caused by two targets. 


=i 





In general "ne = T) aoe Ab 





2 
eet 2 ve 
eee ele 2 2 _ 
= ={[ 5 f S(t) dt] (353-2) 
et 
T +At T5tAt 
Tae Pa sails 
» is f Sejetel eb ji CRG kota 
T574t wasn 


anGazat — An > Q> 


1 Te 
0 ney T,~4t 
ay : -T)74t =< E< TS74t 
SE) = ans T,~At << Vices T +At 
ap ee ee At 
0 ee tet T+At 
+ 
qT, ie 
A 
y= f SiGe nate 
T,74t 
x = Bie! 
a Ja, ta, er | 
oO 1s the relative phase between the two 


returned pulses. 


l is the position of the true 


target. Thus we get an error in the computed range. This 


error is between the two extremes given by the cases when 


the two pulses are combined in phse and out of phase. When 


the two pulses are combined in phase we get a simple 


Beeeression for (3.3-3): 
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AL = ————— - T ee 4) 
2 


The error due to the addition of the second target is 


given by (for the case that they are in phase): 


which is not equal to zero when T, = T.- 


oe CONCLUSION 
In the last two paragraphs we have shown that when 
there 1S an additional target in the resolution cell, the 
Fadar will not track on either of the targets, but on 
their center of gravity. In order to eliminate this 
phenomena, we must change the design of the radar so that 
it will track on the two targets. 
There are two main ways for accomplishing this objec- 
tive which have been developed in this study. 
moe extract information of the location of the target 
by processing the signal after the recevier. This 
leads to the implementation of what we have called 


the probability filter (Chapter V). 


1 ei = ——— (T, -T,) OS Sas, 


(11) Modify the receiver and the filtering process after 


mem Gecelvyes in Order to extract the information of 
the two targets. 


This leads to an antenna modification and a new type of 


modified Kalman Filter (Chapter IV). We will consider later 
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the case of a ship as the true target and chaff as the false 
target. This implies very small Doppler shift between 

the two targets. In the problem of air targets the differ- 
ence in velocity between the target and the chaff is 
distinguishable anal they can be separated easily by a 

moving target indicator (MTI). In summary, the following 
techniques are presented in this dissertation: 

(1) Antenna modification plus a Modified Kalman Filter 
(M.K.F). In this case the location of the targets 
are resolved on a pulse to pulSe basis (See Chapter 
IV-A). Thus it can include the case of air targets 
where doppler shift is not applicable rather than 
be restricted only on surface targets. 

Bee CObability filter. This solution is suggested in 
case of C.S. or L.S. radars. Here the location of 
the targets cannot be resolved on a pulse to pulse 
basis (See chapter II), and many pulses must be used 
in order to solve the problem. Because the dynamic 
of a surface target is slow, we can assume a 
Stationary process for the target's RCS. Thus a 
time average over a small interval can used as a 
good estimate for the mean of a function of the 


Gre 
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V. A MODIFIED ANTENNA AND KALMAN FILTER SOLUTION FOR 

UNRESOLVED TARGETS — FOR A MONOPULSE SYSTEM 

fSweeethe, THEORETLCAL SOLUTION 
fm Introduction | 

The problem of two targets in the same resolution 
cell is discussed in the literature, mostly from the point 
memvtew Of resolution and multipath [Refs. 12,13,14,15,16, 
MmeeeGo, 19, 20]. The principles involved in the multipath 
problem are Similar to the two-target problem, and it is 
proposed that they be adopted for this case. The major 
difference between the two problems 1s that in the multi- 
path problem the amplitude ratio and the relative phase 
between the two return Signals are known. This information 
is missing in the two-target problem. The method used to 
solve the multipath problem is called the "Complex Angle" 
method (CA). In this chapter we resolve in a new way the 
positions of the targets in a single pulse and demonstrate 
how to filter the data. We break the problem into two 
earcs: 

a. Resolve the poSitions of the targets and their 
amplitudes, assuming that there iS no noise. fThis 
gives us the positions of the targets plus estimation 
noise. 

b. Filtering the positions of the targets. Since the 
foregoing yields the position of the target plus 


estimation noise, an additional filter must be used. 


oy 








Furthermore, the estimate of the positions is with- 

out the decision of target vs chaff. Thus this 

filter must select the two positions related to the 

eMart and the target. 
The general assumption is that in the input to the antenna 
there is a white Gaussian noise. Since the first block 
includes a non-linear transformation of the signal (see 
fig. 4.1-1), there is non-GausSian noise at the output. 
ieee output is the input to the filtering block. Our 
filtering procedure is complicated even with the assumption 
of a GauSSian process. Thus, to save a lot of computation, 
the assumption is made (Similar to the usual monopulse analysis), 
that all the variables are Gaussian. In regular monopulse 
systems the measurement process requires Reece te ataien 
A) {see II.B). The variables A and } are assumed to be 
Gaussian. Hence the meaSured error is not a Gaussian varia- 
ble although in most of the regular monopulse radar, A) is 
assumed to be GausSSian variable. The same type of assumption 
is made in this application. 

ivewerockiealagtram 1S given im Fig. 4.1-l1. We 
start with the first block, i.e., the resolving procedure. 
feescoolving Targets with Monopulse System 
This paragraph introduces a new technique for 

resolving two targets, based on the ‘eoma les: angle" method. 
The uniqueness of this solution is that we use only a four 


element antenna, which is the regular monopulse. However 
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we require four output channels from the antenna rather 
than three channels as in regular monopulse. The resolved 
positions, Similar to regular monopulse, are made in one 
pulse. 

For monopulse system and two targets, the output 


Seeeach antenna is: 


e. = AG, (x, /Y,) + A,6; (%5/Y9) ie ee (ae J.) 


meen vector form: 


e = A,G(X,) + A,G(X,) 


where 
Fela ’ 
G, [wclie vo. tage pattern Of the 1 recelved 
beam 
ee coordinates of target a: (X7¥ 0) el 
A, - effective complex voltage due to target < 
(normalized to the range), and is given by 
Ns Ko G[X ,lexptj lw, t ~ opal? (4.1-2) 
eee So the doppler frequency of target a 
8, 7 phase associated with target a 
2 Jol 
oz The a target voltage cross-section 
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G - one-way voltage pattern of the transmitted 


beam 


K - normalized factor. 


Note: For implementation, it is much simpler to assume 

a) = O and then solve the equations. This is reasonable 

if we deal with surface targets, and so ei 0, and this 
assumption would simplify our decision procedure. However, 


the reason that we do not assume Yy = 0 1s because 


oy is a random variable with almost zero mean but with 
variance not equal to zero. For example, a ship is tilted 


by the waves and the c.g. of the ship is changing with 


time. Thus, the assumption that Yy = 0 can cause an 
error. Also, the assumption Y = 0 would lead to a special 
case solution which is not desirable. 

There are four unknowns for each target: 

= The positions - xy = (XX) 

- The amplitudes - Dees (Re Avr Im AY) 

Because eq. 4.1-1 is a complex equation and we have 
four antenna elements, we have a total of eight equations. 
Meter SOlving (4.1-1), one can use (4.1-2) to solve for Oy 
and the angle (t° Wa, + a): By uSing at least two different 
Samples, we can also separate Way from an and have the 
doppler information about the target. The main approach of 


this section is the solution to (4.1-1). Because this 


equation is quite complicated, the main idea is: 
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mee solve it by approximation, using Taylor series [12] 
as is usually done in regular monopulse analysis. 
(2) Solve the non-linear equation with -a starting point 
given by the first step. 
The Taylor series solution for this case follows 
with the solution given by (4.1-18,19). 
Our assumption is that the beam is composed 
M@eecactors, each related to a particular angular 


eeemension, i.e.: 


G. (X) = KX, (x) °¥, Cy), ee ed 2, 4 eae 3)) 


We solve the equations for the unknowns by expanding the 
gains, X, and Y,; by Taylor series up to first order about 


the boresight line, i.e.: 


X, (X) = eel + 4, X) + 0 (x*) 
(4.1-4) 
Y.(y) = vg (1 + By + O(y*), 4 = 1,2,3,4 
where 
g - is the boresight gain 
eS is the slope of the beam near the boresight 
[pcm OWne cs ciuect 1 On 
Bs - is the slope of the beam near the boresight 


lbmMewbOrey Glrectlon. 
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Meom (4.1-3) and (4.1-4) we get: 


le oe, o 


Equation 4. 1-1 can be approximated by: 


Be SGP (Le Cha f(s Soa a (4.1-6) 


+ GA, (l+a,X,+8,Y, +4,8,X%,7)) aioe eh 
Those are four complex N.L. equations. Let us now 


"linearize" the equations by substituting new variables: 


JE aL 2 

4. = ALA aT A,X. a a a 
Z. = ALT, + A5Y. 

Zi = ALA 24 oF AjXjt, 


In this case the linear approximation of (4.1-6) is: 


Pe 2 eee oe 2,)-G , i = 1,2,3,4 (4.1-8) 


67 





earn Short form: 


e = Ba. ©. (4.1-9) 


where: 


| 


| 


g ; ° ; ; (4.1-Sa) 


We can assume that Bt Gxrsts ana sO,cne SOlutIOn of 2 is 


given by: 
Z = B Gc (2B ks) 


From the knowledge of Z we have to find the location of 
the targets (X), and the amplitudes. By equation 4.1-7 it 


~~ 


is easy to determine the amplitudes, given the locations: 
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a 5 LO a ae a 
al X,Y, = XoY) X57, = X,Y, 
(411) 
A : Son. ete 42%, : ~4, 4,7, ar 44 
Z, X,Y, = X54 X57. = XY 
: Toe 
Let's now define new variables (a = [a,,a57a,7a,]): 
3 A X Xo (Y, = Y,) 
1 X51. = ALY 
ae ees 
Z XoYo = X,Y, 
Cae L=i2) 
. A Y,¥5(%, X,) 
S X51, = XY, 
a A ae 
4 X5Y. = XY) 
This gives us a new linear form of equations 4.1-11: 
4,44 aP Z4a5 = 4 
(4 ee) 
4,43 at 442, = 43 


The last two complex equations (4.1-13) are four equations, 
and we have to write Z;as it is the sum of real and 


imaginary parts. 
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este, f= 1, ..., 4 (4.1-14) 


Let's define 


ic 


tH 
ic 
-—__ 
coe 4 
fx 


where: 


We now can return to eq. 4.1-13 by defining: 


ic 
be 
Re a 
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Solving this equation gives: 


> 


Assuming H 


Z 
solve X.- Murserwe NOtlce that: 

a ee 

=A ee 

oe 

aes a 

a) 

By SE 

- Se 

2 iy aT 


can 


exists, then a is determined. 


in) easy 
(Beet RT) 
(-R_I., + RI.) 
avi lve a Leakey 
ae er TS) 
(-R,1) + R15) 


Now we have to 


Coes ) 





For convenience we define a new variable: 


U = HX 
be We 
7 
Y 
e 2+ (4.1-18) 
oD 
Ly] 
| 
0, | 0 
H = ————4 
ad 
| 
-1 i 
A 
D, = 
i i 
The inverse transform is given by: 
x = 2H U OT ig 
me eerie es ; 


memving (4.1-17,18) for U gives, by eliminating xX: 


UG 





2 
a, Se 
a,a, - a2a, + 4 
aA 
AU Pies C210) 
a 
S Se yh — 
2 
4 - a1, -- Aza. 
25 


when S 1s a variable which can be +l or -l. In determining 
U one solves for U, and U, first. Using these values one 
een find U, and U,- Piomioe Tl Lustratked Im Fig. 4.1-2. 

The variable S occurs because we have a nonlinear 
equation. It is easy to see that changing S from +l to 
-1l will cause the whole solution to change indices. (lI.e. 


given e 4.1-10 yields Z, R, I and b 4.1-16 yields a, 


De 
fed Gives U, 4.1-19 gives X, 4.1-11 gives Ay and A, 
see figure 4.1-2): 


When S = lL 


3 





*sjosie] ayi 
Jo suotytsod ayy pue sopnazty{due 





ay. Bulutusaiyap Jo weirsetp YoolG “72-Po pr aansty 
Y <j ae 
404 7 
BuLALOS 
NG /\ 

(Qs | peubls 
~ voIUadNasay 

ERE) Ray 1G 

wu< i 

X 


| = 

= € Se 

Ta %o Lat} 

ji AO J ae ee Oyauedas 
; Z 3 


aes) 
eS EE, 





aseud 


BULA 
WAopJSueul ii Ues. 


= 
[ 


n 


CT 






Sas|[nd sapeubis 


paztuouydudrs sSeuua UY 





Xo 
Bos 
ae = é 
al 
al 
PAL 
fbx, = 


: 
| 
aed 


From which it is clear that in both cases we have the solu- 
tion of the location of the two targets. For the solution 


(Sy) assume an arbitrary that S = +1. 


“) 
il 
iL 
A i A 
ia 4 — a a — ~1 
Meso lution (2) 

for S=l i Ko 
% 

; 
i As 


Then we have a unique solution for the positions of the 

targets but we can't match the position of the target to 

the measurement. There are two possibilities for target #l 

position (and the same for target #2), i.e.: the position 
il eP 


of target #1 can be (.) or (."). 
Yi ¥ 2 





The block diagram of solving the positions and the 
amplitudes of the two targets is given in Fig. 4.1-3. 
Everything is straightforward in the block diagram except 
the amplitude vector for the range, which will be explained 
in the following paragraph. 

Solution for Range 

AS was mentioned in Sec. II.B, the error in range 


1S given by: 


Vie Vo 
= x —_ 
AR K pay (4.1-2))) 


where: 


ces - 1s some measure of the amplitude of the 
X gate (early or late) 
Kr pee Ol) oS Gallia. 


In our caSe we cannot take the amplitude in each range 
gate and substitute into the formula because at each gate 
we have the sum of two amplitudes. Thus we have to use the 
technique developed for resolving the angle. The simplest 
way is to use the procedure for solving the amplitudes for 
each target, but now we have to consider two amplitudes for 
each target (for the early and late gates). The procedure 
is described in the block diagram 4.1-3. 

(1) For solving the regular amplitude we use 
formula 4.1-11. We need the delay because we need to wait 


for the solution for X and Y. 
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(2) We take (at least) two samples of 2% in the 
late gate and in the early gate (2,2) and we use the 
same formula for obtaining the amplitudes at each gate, 


where: 


For the purpose of tracking in range we do not need the 
complex signals. Thus we give the output only for the 
absolute value of each amplitude. 

Note that here we cannot use integration as sug- 
gested in Sec. II.B because we have to take a sample of 
the amplitude. To perform a process like integration we 
have to take many samples. Our process involves a phase 
detector plus a computation for each sample, and this limits 
the samples that we can take. In our case in order to save 
the amount of computations, only two samples are taken. 
Furthermore, we can save even more calculations by assuming 
the range of the target and the chaff to be the same. 

After we have the approximate solution for the 
unknowns, one can substitute these values in eq. 4.1-l1 as 
a first approximation and iterate the N.L. equations for 


amore accurate solution where we can assume that: 
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sin é, (u - aa 


= a —— zh 
U, (u) = K ooh = eee 2 ae eo oo. 22 ) 
iL oi 
with: 

Dees tS (XS Ory. 

= ad 1 
u = as the coordinate x or y 

ba K, - parameters of antenna 1 
Jaa - the offset of antenna 1 


Let uS now observe that this method gives us exactly the 
same solution for the single target case. In this case, 


A, = 0 and so (4.1-7) reduced to: 


1 1 
ba ES 
gee. 
Be SG 
i ea | 
And so 
, ee 
1 By Be 
(4.1-24) 
1 By ze 


We) 








Now, let's assume the same beam shape for all the antennas 
ie.e., jos | = 18. | =a). Then the matrix B becomes 


(refer to Fig. 2.2-6 and equation (4.1-9a): 








i or} —O ara 
si = O — 
iL —O -a a | 
1 a O aa 
a 
Betvying for Bee 
2 2 2 Z 
or} O or} —O 
O O —O —aO 
pl = = Pa DIG | 
+90. —aO Qs —o O 
-1 an il: -1 


Now we can solve for Xi Yy (Alea) 2 


(2, E e,) 2 (e, + e,) 


il L x 
x = = = = Ca l=2 7) 
A. OL er ao e, + e4 + en eo } 
RR RE? ek el ee 
aL OL e1 ~ e, ~- e4 ~ en Go 


which is the well known solution for regular monopulse 
Systems. Thus, from this result we can say that regular 


monopulse system iS a special case of our development. 
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Pee ine Filtering Procedure 


In the last paragraph we investigated a way to 
resolve the unresolved targets where we have a monopulse 
system with phase detector. The development was based on 
the assumption that we have no noise. We could work this 
problem as was suggested by [20] but this leads to a lot 
of computations, and the estimator functions on a pulse to 
pulse basis ({20] does not use the prior information). fThus, 
after our estimator we have to use a filter. Furthermore, 
we have also to make a decision about target vs chaff. 

For simplicity we assume that the outputs n (See 
Fig. 4.1-1) are Gaussian variables, so we must solve two 
problems: 

(1) Selection - 1.e., at each pulse we have to 
@eeeinguish the origin of the measured aae a, which means 
that each measurement gives us two sets of data, target 
and chaff, but we cannot distinguish between the two. 

(2) Filtering - After we make the decision, our 
data iS not noise free. Furthermore, the decision may 
be wrong. For this purpose, we have to filter the data 
Output from the selection block. In the previous analysis 
we assumed that measurement noise does not exist. In 


reality we should have started with the equations: 


23 in A,G,; (x,) + A,G; (x,) Be n, i= 3p Ol tes (4 1= 28) 


81 





where n, 1s a complex white Gaussian noise: 


ee es 128 ee (4.1-29) 


Tau, tl ee =—ewhihte Gaussian noLse with zero 
mean, each of them with variance 
O°, assuming no correlation 


between no and oye 


However we used (4.1-1) so that after our processing, we 
Can assume that a GauSSian noise is introduced. The 


computed information needed for Gaussian variances are: 


eis) The means 
Heal: ) The varlances 
q 
ey. ) The correlation between the noises. 


That is, we have to include, for the random process n: 
u = Et{n} (4.1-30) 
oa COVAR{n} (4.1-31) 


Which 1s straightforward. 

Once we have mn output of the resolved block positions with 
the approximated statistics of the noise, we can make our 
decision and filter using the information about the ampli- 
tudes of the signals and the dynamics of the targets. 

The decision and filtering process is described in the 


block dlagram given in Fig. (4.1-4). 
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Our information comes from the monopulse antennas, 
and we compute the positions and the amplitudes of the 
targets in the resolving block. The outputs from the 


ae) (Ay 


resolving block are n balou which include the posi- 


tions and the amplitudes of targets number one and two, 
and we have for the filtering procedure to decide if ae 
belongs to the target or the chaff. The information needed 
to decide is the predicted position of the targets and 
Siemcovariance matrix of the error of our prediction. 

After the selection block, we get the information 
needed for the filtering, i.e., the amplitudes and the 
positions of the target and the chaff. The output of the 
Filtering block is the estimated position of the target 
and may be, if needed, the position of the chaff and the 
velocities. 

There is also an outside output to the filtering 
block, 1.e., during the process we change the position of 
the antenna in order to zero in on the target, Because our 
model refers to the absolute coordinate we have to take 


the change of the relative position antenna-target into 


account. In our case we make the above assumption for 
Eeimplicity: 
1) The process is sequential. This means that we 


calculate the estimator output and the covariance 
matrix from pulse to pulse. We do not have to store 
all the measured information so far. This assump- 


tion comes from the desire for a simple procedure 
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and reduces the amount of storage and computation 
needed for a non-Sequential process. 
ii) Our variables are Gaussian. 
The theory of the selection and filtering is given 
in Appendix A. We give here a short description of the 
procedure. 


at) The dynamic model: 


A{(K+1) = OX(K) + AU(K) + W(K) 
x fT) ox) 
X(K) © | ------ = the amplitudes, (4.1-32) 
a A and the 
cme 
(Cc) positions of 
| X (K) | tne targets, 
~ Xx 
ae 
_ 
g (7) | 0 
| 
7 is, ae 
— | 
| 
0 3 (C) 
| 
aed = the transition matrix of each target 


(1=T or 1=C) 
a 

ee pees. oe 
| iS) 
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11) The meaSure model: 


Me LIS) C X(K) + V (K) 
Gre We eae <)) 
Te (K) 
n(k) & | ------ 
re) 
Pin is Ene vrandom Switching matrix 


(see eq. A-l). 


111) The statistical information about the states and the 

noise 1S given in Appendix A. 

a. The Selection Procedure. All the procedure is 
described in Appendix A. Hence only the block diagram is 
given here (4.1-5). 

Description: We have in this block the above inputs: 
i) The resolved information about the targets - n (K) 
Hei ) The estimate predicted positions and amplitudes 
of the targets - ¥ (K|K-1) 
111) The covariance matrix of the error of the estimate 
prediction - A. 
The output from this block is the assumed information about 
the target. The resolved information of the target is 


transferred into two forms: 
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aio 
i) t, = = which is the original input. 
fa 
n (4) 
11) t, = ---- which is the switched input. 
(1) 
J 


After we form the two possibilities ti, we subtract the 
Peeaicted measure ¥(K[K-1) with the formed t, and we get 
a two error vector r,. We now form the likelihood test, 
and we decide t, or ty according to hy > Lo Oe Lo < dy. 
The output from this block is the assumed information about 
the target. 

bee Filtering. 

After the decision about target vs chaff is 
made, we have to filter the output from the selection block 
(4.1-5). The filtering could be straightforward if we have 
moiecontidence that our selection is correct. In our case 
we are not sure if our selection is correct, and we must 
take this into account. Since we cannot tell if our selec- 
tion was correct or not, we can only take into account the 
probability that the selection waS correct, 1.e., we cannot 
reject or accept the data, we can only "weight" the data 
depending on the computed probability that our data is 


merrect or not. 
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The two possibilities are: 
a) A true data available, i.e., we decide the right 
selection. 
ii) A false data available, i.e., we decide the wrong 
selection. 
We must aSSign a certain "weight" (depending 
on the probability that this data is true) to the data. 
Because all the procedure is descibed in Appendix A, 


Seer tne block diagram is given here, (Fig. 4.1-6). 


Bb. SIMULATION OF THE MODIFIED KALMAN FILTER 
fo Lntroduction 

ime M:K.F. is very important in our system, and 
we must show that the concept of the modified filter works. 
We will use only one component of measurement to show how 
it performs. Because the most distinguishable component 
between the chaff and the target is in the height, the 
Simulation is only for height. To take into account the 
Variance and the bias of the estimation error the probability 
of missing is used aS a criterion. The probability of 
Missing includes both the bias and the variance of the 
estimation error. Furthermore the radar system is assumed 
to be in a missile. 

We have chosen only one component of measuring, 
1le., the height. The target which is assumed to be on the 
Surface ,=has zero net velocity in this direction. The chaff 
has no negative velocity in height. Thus we can define, 


Semeeocaing to (4.1-32): 
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eae (1) 7 . 
X, (K) yf) (K) | 
| t 
=“) 2 |x| a A ashlee | 
Z ey te (0) | 
j 
X., (K) x 6S) ogy, 
x, (K) . = 
a SL 
where: 
X, (kK) = xT) (xy - is the height of the 
target at time K 
X. (K) = xO) (x) - is the height of the 
1 hati. at time K 
X. (K) = Ko - is the velocity of the 
~ chaff at time K 
and 
ios) 0 


According to the above assumptions, we can relate (4.1-32), 
(A-2) and (A-2a) with the proper matrices. The matrices 
are changed according to the simulations performed. 

Note: In the simulation we use the terms (refer to (4.1-32)): 


- "without switching" which means 


L(K) = ceorna link: s 


- “with switching" which means that L(K) 1s 


aSsedertnecdelm. (A—1p) - 


eyaL 





fieene Simulations we compare the M.K.F. "with switching” 
memene R.K.F. “without switching" as a referance (ideal), 
and to the R.K.F. "with switching" to show the improvement 
that we get. 

A Monte Carlo simulation with 200 ensemble members 
is performed. The results are given in comparison to the 
regular Kalman Filter (R.K.F.) (where there is no switching 
M@mehne input). The criterion for missing is as follows: 

- The missile is Simulated aS an inertial mass 

and it is guided to reach the estimated height 


Srernemeadrget, 1.e6.% 


ey, 


Xy, (141) = KX, (2) fp (1-K) x (ales il), Care 2) 
where 
Ky (2) Me ecemmssinle SemOsSltion at time 1 
x (T) ay - the estimation of the height at 
time 1, see (4.2-1) 
K - transition constant of the missile. 


Tt is chosen to be 0.8. 


Ztie Catterton Of MLSSing: 


x'T) (pre) | < o Ueeaae 


[ES - 


where: 
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TTG - time to go (chosen to be 40sec) 


(1) - the height of the target at time 
1, see (4.2-1). 


In our case 8 is chosen to be 0.75 which experi- 
mentally gives a probability of missing of about 


Ono 


We study the performance of the two filters by changing 
three significant parameters: 

- variance of the noise. 

- mean of the initial states. 

- standard deviation of the initial position of 


Mem Cian i je WiC. 


= 2 
aay ENE (O) = 000) 1.2} 


meee Ln the Simulation instead of the formula (A-13) for 


a(K) the following formula was uSed; by mistake: 


a (kK) = Beer LWP A a. 
: 1 + y(K-1) 8 (kK) 
where: 
y (K-1) joe 
tb & pl cx) a7 (KR) bl x) 
SOR one x (RK 1) 
ee “1 | *2 
oy US Sal) ae (IK | egeee 


and §(K) is the same as 8(K) given in (A-13a), and 
A+ (K) is defined in (A-6). This leads to an a(K) 


Slightly higher than the real a(K). Apparently this 


She 





elangewim (kK) will not affect the simulation results 
Significantly. Hence it should not change the 


SmitlarLon Conclusion. 


2. Simulation Results 
a. Influence of the noise on the Performance of 
the Filter 


The following parameters and matrices are chosen: 


Ih 0 0 
o = 0 ue 1 
0 0 IL 


pe 
II 
aa | 
ro) im 
be ro) 
(ae 


ag 0 
M = @) 5000 0) 
0) OEe5 
10 
X(0) = 100 
ack 
a 0 0 
es = 0 a/100 0 
0 0 a/4 
where 
"a" -—- ais a parameter to be changed. "a" 


is changed from 0.1 to 100. 
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Three results are given here: 
meee fine probability of missing as a function of -1010g (a) 

(Fig. 4.2~1). From the Simulation results we see 

mma t : 

() The R.K.F. (without Switching) gives a better 
result than the M.K.F. (without switching), 
but they are close. Thus the MKF approaches 
the ideal. 

(11) The M.K.F. gives significant improvement over 
the R.K.F. where the input to both is switched 
randomly. (The probability of a miss for the 
R.K.F. in this case is almost one). 

(2). The difference in gain between the M.K.F. and the 

Ewes. aS a £UNCELOn Of —-1O log (a) (Fig. 4.2-2). 

The comparison of the gains at TTG is given here. 

AS was expected, the gain of the M.K.F. is lower than 

the R.K.F. and reaches the regular gain as a +0. 

This iS an example where the M.K.F. approaches the 

Beker. im the limit. 


Di. Influence of the Mean of X(0) to the 
Performance Of the Filter 


The height and the velocity of the chaff are 
taken in this case aS parameters. Both are distinguishable 
parameters between the target and the chaff. Thus we 


Study in this simulation how they influence the performances 
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of the filter. The most distinguishable parameter is the 
initial height. The following parameters and matrices 


are taken for this test: 


1 0 0 
og = 0 h uE 
0 0 i 


50 0 0 7 


1.0 0 0 
Q = 0 9-107 0 
0 0 On.25 
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where "a 


1S a parameter to be changed. 


Two results are given here: 


(ey) 


The probability of missing as a function of 
ae (X) (0)) (see Fig. 4.2-3). We see 

from the results that the probability of 
amiss is larger than that of the R.K.F. 
without switching, as can be expected, and 
it reaches the performance of R.K.F. in the 
limit (when at the instant TTG the two 
heights are far apart, a > 100 or a < 10). 
The maximum aieeineinl Lien of missing is 
reached at about a ~ 50 (target and chaff 
aS approximately at the same height at TTG). 
The initial gain as a function of "a" in 
Bac 4 2=4 ie Here we see that in the limit 
the gains are the same as for the R.K.F., 
and reduce to 0.5 when the two heights are 


indistinguishable as expected. 


The following inputs for testing the influence 


of the velocity are now chosen: 
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where: 


iL 0 0 
6 = 0 il iL 
0 0 1 
al 0 0 
Cc = 
0 il 0 
i 0 
R = 
0 i 
ibe 0 0 
Q = 0 S510) ante 
0 0 Be 
25 0 0 
M = 0 25 0 
0 0) ae 
nO 
con = 10 
a 


is a parameter to be changed (the initial 
velocity of the chaff) 


Paz 





eee eect OO 1 


c = 2b 


Note that in this test the mean of the initial height of 

the target and the chaff are the same. So the only dis- 
tinguishable part is the velocity. The target has velocity 
zero, so the most difficult decision is at velocity zero. 

In order to check the case when the target and the chaff 

have the same dynamics, we choose b and c so that for a = 0, 
the chaff remains at almost zero velocity. For this case 
(chaff velocity about zero) we must expect about 0.5 proba- 
bility of missing. As the difference in the two velocities 
increases, we can expect that the probability of missing will 
approach the probability of missing for the R.K.F. On the 
other hand, for the R.K.F. where there is switching in the input, 
the probability of misSing is at the lowest value for 
velocity zero (because both of the targets have the same 
height), and approaches one as the velocity varies from 

mero. The results:of the Simulation are given in Fig. 4.2-5. 
The results agree with the discussion above. 


es Influence of the Standard Deviation of the 
fiteiainotate tothe Performance Of the Filter 


The standard deviation of the initial state of 
the chaff changes the probability of the decision at the 
beginning of the process. 

Thus it is interesting to see how it influences 
the performance of the filter. The following matrices 


and parameters are taken for this simulation: 
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Be 0 Vie 
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KD 100 
mee 
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M = 0 a’ 0 
0 0 ORS: 


Where "a" is the parameter to be changed. 
The difference in probability of missing, as 
a function of "a", between the M.K.F. and the R.K.F. (without 


SWitching) is given in Fig. 4.2-6. We see that in the 
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limit (a > 0), the two performances are quite close, and 
when the standard deviation increases, the performances are 
Baitedpart. iI.e., the probability of missing for RKF is 
better by 0.1 over the MKF for a = 100). The difference in 
gain between the two filters is given in Fig. 4.2-7. [In 
the case Go, (0) 1s most influenced by the parameter and 
we Show this difference in Fig. 4.2-7. As expected, the 
difference between the gains increases as the standard 
deviation increases. 
ad. The Gain as a Function of Time 

To illustrate the behavior of the gain as a 
Function of time, a typical example is given here. We know 
that the gains, contrary to the R.K.F. are data dependent; 


so for different data we get different gains. 


The inputs to the system are: 


HE 0 0 
og = 0 Le L 
0 0 At 
al; 0 0 
ee == 
0 i 0 
Ji 0 
R — 
0 1 


108 


—_ 








Q = 0 9-10 0 
0 0 One 
10 
X(0) = | 100 
~1 
50 0 0 
M = 0 10 0 
0 3) Ons 


Pemexolanation of Fig. 4.2-8 folléws. At K = 1, a is about 
Moye enus the gain is reduced to about .5 of the R.K.F. 
gain. At K = 4, a becomes almost one and the gain first 
increases and becomes larger than the R.K.F. gain, and then 
it decreases below the R.K.F. gain because a is about l 
(a = 0.998). It then reaches a steady state gain. At 
K = 35, a is reduced to a low number (a = 0.57), because of 
the large noise in the input. Hence the gain reduces to 
almost half of the steady state gain. At the next instant 
(K = 36) © becomes again about 1 and the gain is above the 
R.K.F. gain. After that % remains almost 1 and the gain 
decreases again and reaches the steady state value. 

For the M.K.F., the steady state value is 
reached only for the intervals where « is almost one for 
Beveral K's. it can deviate from this value for short 


transients because of spurious noise. 


ee a 
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Bee Conclusion of the M.K.F. Simulations 
Based upon the simulation results we can conclude 
ema t : 
a. The performances and the behavior of the M.K.F. 
is as expected. 
b. The M.K.F. is applicable to the chaff-target 
problem and handles the ambiguity problem in 


the monopulse solution (see Chapter IV.A). 
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econ! TON HOR UNRESOLVED TARGETS FOR C.S. OR L.S. 
meme RODABLLITY FILTER 
Ae THEORETICAL SOLUTION 
In case of L.S. or C.S. it can be assumed that a phase 
detector is not available. Thus we get only amplitude 
from the receiver. The approximated p.d.f. of the return 
Signal exists in the literature [8], [9]. Our purpose is 
to approximate the amplitude of the target or the chaff as 
received at each antenna so that the position of the target 
can be calculated. In case of C.S. we can separate the 
Signals received during one conical scan period as equiva- 
lent to n Separate antennas where n can be as large as the 
number of pulses in one rotation. 
The information of the amplitude of the target at each 
antenna is sufficient to determine the position of the 
target with respect to the antenna (see Chapter II.B). 
The probability filter proposed here seems to be the 
MOSt practical solution to our problem because: 
(1) It gives a good estimation for parameters 
(Signal amplitudes in this application). 

(2) It is easy to implement. 

(3) The probability filter does not depend upon an 
exact knowledge of the p.d.f. of the signals, 
and an approximation is sufficient. This filter 
can be implemented with a premeasured reference 


signal which is used with an approximate p.d.f. 


aIEZ 





to obtain the solution. This is an advantage over 
a likelihood filter which cannot be implemented here 
because the p.d-£. of the Signals is not know exactly. 
The theory for the probability filter is developed in 
Appendix D. The application of this general approach follows 
with simulated examples to demonstrate its effectiveness. 
Because we restrict ourselves to surface targets we can 
assume a Stationary process for the amplitudes measured at 
the receivers. 
Thus we can use (D-12b) as an approximation to the time 
average (l1.e. with reasonable "T" compared to the missile's 


fmrant time). 


Be PROBABILITY FILTER SIMULATION 
mor cntroduction 
In this section we give some results using the proba- 
bility filter concept. The difference between the likelihood 
filter and the probability filter is that for the likelihood 
filter we need more information about the process, 1.e., we 


eed to know the joint p.d.f. (f oe while in 


y(SprXorees 


~~ 


the proposed probability filter we have to know only the 


Memoonal p.d.f£. (£.. (x)) or its approximations. The relation- 


Xx 
il 
Ship between the samples involved the expected value of a 


function of the samples, and for this the sample average 


is used (assuming an ergodic process), l.e.: 
B{f(z)} = £(2) = = J £(2(4)) (5.2) 
a 
where the symbol "~" means time average. 


ae 





In order to estimate a parameter vector 0 we develop 
the theory of the probability filter given in Appendix D. 
The main idea of the probability filter is to minimize the 
difference, in some sense, between the conditional p.d.f. 
£[2| 6} (where 8 is the stimate of o), and the measured p.d.f. 


of the process, Z, which contains the information of 64, 


£(Z) = £(Z|@). We choose the norm to be (D-11): 
r= + fte,(z|6) - £,(2)17 az > min 
2 BN Z 


am OO 


This leads to a necessary condition (D-12): 


9£, (Z| 8) 


0 = = i sg ptt! ?)az a 


cee (2) i 
gemrcmmez 


Q? 
FH 


Q> 
¢ D> 


Our task is to solve this equation for 4. 


Summary of notations: 


: = estimation of 6 

o = true value 

8, — 3th iteration in the process of finding °. 
8 (K) = : at the end of ries data block. 


From D-12 we want to determine 98 = 6 so that gI/ 96. = 0 
fer all i. In order to accomplish this the following procedure 


is used: (D-11) and (D-12) can be rewritten as 


joe) 


ee’ hearer _ 
be J eae) aya ale (5.2-l1a) 
oof, (ZY) ogf,(Z|Y) 

iy J —“sy—F, (Z| ¥) dz - J a¥ £7, (2)d2 = Eq (¥.8) 
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The experimental value of dI/93Y is given by 


Bea Y) 288 (Aes) | 


N 
oI - A ie eae Z A 
x ee! JY £,(Z|¥)dz n.2 oY ~ £, (7,8). 
(52 2=15) 


Note that £,(Y,9) Goes not depend explicitly on 6, although 


8 of course influences the Z(1). 


nw 


We are seeking @ such that £ (6,9) 


0. Usually we 


solve this iteratively; ae is guessed, and then 0441 is 
obtained from eal = os - SOs 7 8) Pe the ge 0 Whe Sore 1S 

an arbitrary number. If the sequence 0. converges to a, 

then clearly eae 0) = 0. If ¢ is very small, the convergence 


is slow. On the other hand, if ¢ is large we might get 

oscillations. Thus we must use for this method an optimal 
weeeich might be a function of f.- 
solve for 6 is to use Eqs. (5.2-la) and (5.2-lb) that can 


A different way to 


be iterated so that 


cas nw nw 


Beocme, os) = £..(9. 


cee oy 5 aes ST eae seer ad) 


where a is the pte ape so <itat Lon fOr taba. (5.2-l1c) 1s 


solved successively for 8 


Gaal! Rene glk shefei oer ogee 


“~ 


£3 (8, 78) = 0, and we conclude that y = 0 = 
After we iterate (5.2-lc) for the first N samples of Z(i) 
we get 6(1). Then we take a new sequence of N samples and 


(5.2-lc) is iterated again to produce @(2). Thus after the 


probability filter we get a sequence of the estimated parameters: 


iy egi2), .-.. 8(R), «--) 


ee 





iomactually 0 (K) changes between blocks of N samples, we may 
want to smooth the estimates 6 (K) with a Kalman filter. In 
the discussion which follows the following terms are used: 
- "Open loop" means the eStimation of 8 (K) with a 
Single block of N samples Z(i). 
- "Closed loop" means the estimation of 8 (K) with 
the “open loop" followed by a Kalman filter. 
As we know (see Appendix F), the outputs from 
the Kalman filter are 
8 (K |) - the update estimation of 9 (K) 
6 (K|K-1) - the predicted estimation of 6 (K) 
The concept of the "closed loop" is discussed and 
Simulated for the estimation of one parameter (parts 2 and 3). 
2. Open Loop Performance ~- One Parameter Estimation 
In this section we give the open loop output for 
estimating one parameter. We choose the exponential 


Meet ibution. The p.d.f. is given by: 


emca|ie= che 9 = 2 6 (5.2-2) 


9 is the parameter to be estimated. The equation for the 


probability filter is as follows: 


of ZY) _ zi 
oa ee oe (5.4-2a) 
oY 
ogf_ (Z| Y) ae ° = 
oe Ca 4) az = if Ye ZY az - if ¥“Ze aX ag = 022) 
Y Z 
0 0 0 
(ED Bs) 


Thus we can rewrite equation (5.2-l1b) for our case by 


Meeeteucing (5.2-2a) and (5.2-2b) into (5.2-1b). 





N 
A al « a 
£,(¥,9) = OY! — me “FT ee eee tee) (S23) 


The first term in (5.2-la) is given by: 


Of, (Z| Y) =e aa 
of, (zlayaz = f se (9%) 4az- f aye (OF¥) Zaz 


00 0 0 
_ 2 
= {6/7 (6+Y) | 2-4) 
Meeercnat (5.2-2b) 1S a special case of (5.2-4). From 


equations (5.2-265) and (5.2-4), (5.2-1a) can be written 


for our case as 


te 


£,(¥,9) 22) = 6.25 Says) (8.225) 


New we iterate our solution according to (5.2-l1c) using 


Meee 5. 2-3) to (5.2-5). This leads to the solution: 


cee = G4, fas/(1-a.)] (5.2-6) 


where: 





w~ 


jEtlg Bice) Gaggia S| 





N 
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i=l 


nly 





Note: If the argument under the root turns negative then 
one may take a larger N. 
To show that this concept works, we now test the "open 
loop" criterion. For this simulation 6 = 1 is taken. We 
perform three simulations: 
aeesting —oLl/ay 
The purpose of this simulation is to show that 
the experimental 3I/dY contains the information 
required for the iteration. 
—eeestimating @ in an open loop. In this simulation 


ws 


o. 1S iterated and we try to find: 


(1) The convergence of the iteration by 
(5... 2-6) 2 
eae) The number of iterations needed for 


estimating °@. 
(111) The variance and the bias error of 
é. 
peeeStimating of G@ in a closed loop. 
The open loop estimate of 8 acts like a meaSuring 
M@evice for &. The R.K.F. is used here as a filter for 
this data and we test the performances of the overall 
System. The design of the R.K.F. is influenced by 
the variance of the error and the bias in the output 
Beem the open loop probability filter, 1.e., we assume 


that after the probability filter we have a measured 


Signal with noise. 


| ies) 








otk = O90 (K) + V(K) or 7 cl) 


where 


8(K) is the true value, and 
V(K) is the measurement noise at inStant kK. 


In this case the statistics of V(K) are a function of 


@(K), and this function can be evaluted approximately by 


as 


assuming that 9(K) = 8(K). We are interested in the 
Variance of V(K) which is denoted by R(K). Thus, 
eo VARIV (cK) i=] R(O(K)) = RI6(K)]) (5.2-7b) 


The process @(K) is assumed to be described by a linear 
difference equation similar to (F-2) and the statistics 


mermebar to (F-3,4,5). I.e., 


6(K+t1) = 60(K) + W(K) 
(S56 2=-7¢e) 


@(K) = 98(K) + V(K) 


where, in (F-2), §(K) replaces X(K), and @(K) replaces Y(K). 


The R.K.F. can be implemnted only when 


ie 


R(K) = R[98(K)} RE (K)) 


M, Q, and 6(0) are given. 
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a. Testing the Experimental 3I/dY 
The purpose of this simulation is to show that 
@eecrecOntains sufficient information to close the loop by 
Beeration. 

In this case we take the following parameters: 

(pate = 1.0 

(2) Sample size - 1000 

(3) Number of ensemble members = 1. The 

results are given in two figures: 

Mo wicemoe col ama shig. 5.2-2 which give the 

Tresules of df/oy ftEor 9/8 = O25 £Or-5%.0 2 
memes Significant to note that: 

(cm) At 5/6 = 1 the error between the 
experimental gradient and the analytic 
gradient is almost zero. 

(Gib) The values of the results are very close 
to the analytic values for all - 

fi) — tae curve £or [sel approaches the 
asymptotic line 0.25 which is the 
correct analytic value. 

The results show that the concept of the 

PeOopability filter works. 
Pe EStrmataon of 2 in an Open Loop 

From the first simulation we have seen that 

miy oY contains information about 8. The purpose of this 
experiment is to study the bias error and the variance of 


estimation of 86 for two cases: 


a0 
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- when _ = §@ we will see the influence of 
the number of samples on estimating 8, by 
one iteration (84). 
- when _ is a parameter, we will see that 
the open loop estimator always converges 
toward the true value of 6. 
In both of these cases the number of ensemble members is 
500. The results are given in Figures 5.2-3 and 5.2-4 
Semen Tables 5.2-1 and 5.2-2. 
Refer to the first simulation, i.e. one changes 
the number of samples, N, of (5.2-3). According to the 
figures it seems that the estimator behaves similar to the 
Mikelinood estimator, i.e.: 
- The estimator seems asymptotically unbiased 
(see Fig. 5.2-3). 

- The variance of the error reduces inversely 
with N (see Fig. 5.2-4). 

Refer to the gnd Simulation, i.e., we change 


“A 


oa aS a parameter, and look one step forward to see the 


meror Of oy: We define a new variable: 


wn nw 


a(j) & 05/8, ei ihe eae (sees 


The number of samples, N, for this simulation = 1000. 
Refer to Table 5.2-l1. From the table we see 


that the open loop estimator gives a very good estimate of 


Zs 
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THE ERROR AND VARIANCE OF ESTIMATING 


6. WHEN ae IS A PARAMETER 
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THE INFLUENCE OF THE ABSOLUTE VALUE 
TO THE RELATIVE ERROR AND VARIANCE 


Absolute Value Relative Value 
Q E{8- 6,} var{8,} ai oat)? Var (8,/87) 
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the parameter, even if we are far away from the true 
parameter value initially. Thus it seems that a few itera- 
tions are good enough to estimate the parameter since the 
Single iteration used here gives such accurate results. We 
can use two criteria for determining the number cf iterations: 
a fixed number of iterations or a flexible number of itera- 


tions. For example, a stopping criterion might be 
\150 = By 06) (5.2-9) 


where 





ee lo) 


In this example, it seems that a fixed number of iterations 
are enough. A simulation is now performed to prove this. 
We start with - = 6 and calculate a for a 
range of values of 9. The results are given in Table 5.2-2 
and indicated that the relative variance of 94/8" is 
independnet of the value of 8. The variance is about 


= 3 


2-10 for all the tested values. And the relative error 


is a very low number and it is of the order of om Thus 
we can conclude that our results are generally independent 
of the values of 6. 

To prove that a fixed number of iterations are 


enough (two iterations give a very good estimate), we make 


EAS 





Beweral runs with different initial conditions. The number of 
iterations is 10 and N=1000. The results are given in Table 5.2-3. 
From the table we see that there is no Significant improve- 
ment after two iterations, even if our first value is far 
away from the real one (a (0) = 100 or a (0) =F 0. 0s). Tnus 
for the estimator we use two iterations. From the table 


we see that we get a bias error on the order of 0.23%. 


The relative variance is given by: 


VAR[a(2)] 2 0.002 
Thus we can say that 
a E. 2 
VARI[6.] =— 0.002 6 
and can be approximated by: 
w~ ea n~ 2 
VAR[6.] = 0.002 05 (S221) 


This equation will be used for the noise in the closed 
loop simulation which follows. 
Maumpolosead HOOD Perrormances = One Parameter 
After we have examined the open loop performance 
of the filter, we can close the loop with a Kalman filter 
and estimate the parameter as a function of time. The 
Semceptual block diagram is given in Fig. 5.2-5. The con- 


cept of the filter is very simple. The open loop estimate 
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of 8, acts as the input to a R.K.F. The difference equation 
of 6 1s given by (5.2-7c). The random measurement noise is 
data dependent according to the simulation in the last para- 
graph (5.2-11). So the variance of the estimation noise, 
R(K), 1S computed on line. Hence the Kalman gain must be 
computed on line. 

Explanation to the figure: We use the open loop 
Seeimator for 8.5 aS a measurement of plus noise. The 
initial guess for 6, Dn is the predicted state from R.K.F. 
(i.e. 8, = 2 (Geet). We compute on line the variance 


according to the assumed R(K) uSing equation (5.2-11): 


Z 


Re 5 G@.oae 3, 


(oe 12) 


Nn 


The R.K.F. receives the input @ the computed R(K), and we 


>! 
then estimate the state 8 (K|K) and the predicted state 8(K|K-1). 
In summary; N = 1000 samples of the data are taken, two 
iterations are used to estimate 8(K), so that the Kalman filter 
updates 8 (K) iipuovaol li eyebllter in Eee bEoOckS vor 
1000 samples. The random generators which provide W(K) and 
@(0) are GausSian (see 5.2-7c). Since 8(K) is positive the 
above procedure is made in order to avoid this problem: 

b aL ge is > (0) 

0) _ 

SneGie 5 aL 3e ley 0 
where b is the output from the generator which produces ¢@(0). 
And: 

o8(K) + W(K) if it became positive 


8@(K+l) = 
@(K)/5 elsewhere. 
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The first simulation that we perform is a case where we 


have a process which has the following parameters: 


¢ = 1.0 

On = 0.51 

Oo (O)me=e* 1.0 
VAR{6(0)} = 1.0 


For the Monte Carlo simulation 250 members of the 
ensemble are taken. The performances come out as expected 
with: 

- Variance of the error over the ensemble is of 

the order of 4 - 10°, for all K. 
- Bias error over the ensemble is of the order of 
pee oe) for all K- 
The second simulation is identical to the first one but 
with different parameters. In this simulation we want to 
examine two things: 
- The estimation error bias and variance as a 
PUncEeLOnmot Ke 

- The theoretical variance of the random measure noise 
which is computed according to the experiment formula 
(5.2-12) by comparing the calculated variance in the 
estimation uSing (5.2-12) and (F-10) with the measured 
variance of the estimation of an ensemble of 250 


experiments. 


The input parameters are: 


6 = 0.9 


Q = 0.5 
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femat: 


4. 


6(0) = 5.0 


VAR{@(0)} = 2.0 


Ensemble members = 250 


Mae wBesults (see Fig. 5.2-6 and Fig. 5.2-7), show 


mae Variance Of Ene estimation error is on the 


erder of 2° 10°. 


The bias error is on the crder of owe? 

The theoretical variance of the estimated error 
1s very close to the simulated variance. Thus 
our experimental formula (5.2-12) is accurate. 
We know that our error is biased but it 


is one order less than the standard deviation of 


the error. Hence the bias error is negligible. 


Estimation of Several Parameters by the 


Paceaomihey Eh liter 


iimeroduect 1 on 


imecne Lastmewo SectlLoOns we estimated one 


parameter in an open loop and a closed loop simulation. LF 


Sse section the open loop probability filter performance is 


tested for more than one parameter. For multiple parameters 


it appears that we often cannot find a close form formula 


for iteration procedure for the parameter values, some have 


to solve for them numerically. Performances of the probability 


filter is compared to the likelihood filter. 
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The bias Error. 





b. Mathematical Equations 
AS an example we take a Gaussian p.d.f. with 


two unknown parameters 


H - the mean 
2A 1 
ae 
O 
where: 
2 , ; 
O Ml omene Vethance Of Ehe process. 


The p.d.f. of this process is given by 





£(Z/o,u) = —% exp{- Fa°(2-u) 7} Gay 
V2T 
Let's define the parameters: 
oT ; u 
6 = = 
o. o} 
Sar 4 ai Tl 
a = = (5.2=14) 
eae ) Y 


Z 


Define a new variable related to the real 


parameter (a and u) and the estimated parameters (95,95): 


Whe 





Dei. 2 Z 
B 06 aye! 
A eh eee 2 
= 18, 34 nel ei | 
B 
2 ms ~ 
2 [e.% 6.7 + a? 7] (5.2-15) 
2 iL 
A 2 Zz: 2 
= ©65 Sage: 
9 D 
A AON al OF Ue! STA, aay: 
= J] ace J 5 Vee a 2v5,9, oy 05 
ey, 
8 A a 
_ _ 2 = 2 Ss 2 


The equation of the probability filters follows 


Beeerding to (5.2-1la,b,c) 


maus from 





(522-15) 





aa 


lw 2 2 
(Z-Y,) expt- x¥,  (4-¥,) } 


2 y i e 2 
Lh > By (a, ) SESE Geen 


fon this case: 


By 








(ge Sy ie ; 8£7(Z(i) |X) +¥5) 
oY y E N Han ty) 
ay) = = 
1 N df, (2(i) |Y,,¥,) 
aI a 
[sy7] = Set dY 5 
De 4/1 


(S216) 
The theoretical oT / oy. comes out to be, 


meeerding to (5.2-la): 








y 3 

dT Z 
lay— = 2 exp[- Fo ite 5 

alee V2T 

fn (XY, 9) = = 
¢ 

-) al 
oo — exp [= xb] °d 

25 4/n V27 + B 

(oa2- ly) 
We cannot find an iterative formula relating See Ma ay 


EO O and Me directly. However, we still need to be able 
to solve the relationship 
£, (8,8) = £(8, 9) 
—>£, (6,8) = 0 (5.2-18) 
g= 8 


Assume that we are close to the real parameters and then: 
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Y = u 

b = QO 

oo 

B /2 
; os 
oi i) 
20 


We estimate the parameters in the experimental 
equation as suggested in Appendix D by uSing a reference 
Signal (see Fig. D-5). 

Our system measures [oI/oy.],- The analytical functions 
[ot/oy],, are given by equation (5.2-17). Using the form of 


Meeez-ic), (5.2-18) leads to the approximate equations 











me sel). 24 
J 
(5.2-20) 
we =, 4.2 31 
Osu = a (le a ls 4/7) 


We see that in the limit, i.e., when [ot/ov.1, = 0, we 
get the estimation of the parameters. 
It can be seen from (5.2-17) that the sign of 


the error in uw iS independent of Os and it is dependent 


Only on the error (yu - u) itself. Consequently, the 
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Nn 


convergence is better if ae is fixed and 7 is changed until 
bot/ous zoe Loen M4 and 7 are changed until [dI/oY,] 
and [ot/sy.],, approach zero. In this case it appears that 
we can not use a fixed number of iterations, we have to 
iterate until [ot /oviJ, > Q. 
c. Open Loop Performances 

In the last paragraph it has been shown that it 
is sufficient to treat the open loop probability filter as a 
measurement device and add a Kalman filter to it with measure- 
ment noise R(K) dependent on the data. We perform two experi- 
ments, in which we use the probability filter with a reference 
Signal (see fig. (5.2-9,)). 

The first experiment is to show that the performance 
of the filter is independent of the mean, 9, = un. The second 


I 


experiment tests the influence of the variance, 0,7 on the 


performance of the filter. 
In each of these cases, our ensemble number is 50, 
and N = 1000. For reference we include the performance of 
a likelihood filter for competition. The results are given 
Meerables (5.2-4) and (5.2-5). 
The first experiment (see Table 5.2-4) shows clearly 
that the error in the variance of an and 85 is essentially 
We hold 6 


independent of 8 Pe vecee tt -OUG Case 95 = lO) 


i 2 
and changed oy from -50 to +50. The results for each value 


ohm O45 are very close. 


The parameter 9, is chosen to be zero and 35 is 


i 
meemged from 0.01 to 100. From Table 5.2-5 we see that the 


variances can be expressed empirically. 
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Peak Z(i1) i=l1,.-..,N 









Pocieex(ijt=i,...,N 


Compute £ (94) 










“~ 


Solve for oe aie 








“nw 


£,(0.) = Ep(O5-8544) 





“~ 


Fig. 5.2-8, Block Diagram of estimating 6. 
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aa =e 
2 
(BB) 
A eae 5 2 
CARE) = thas) 3 ale eae 


If we compare the results of the probability filter to the 
likelihood filter (Table 5.2-4,5) we see that the performances 
are close, but the likelihood filter gives better results. 

The experimental results, which are given by 
formula (5.2-21), must be used when closing the loop with 
a Kalman filter as discussed in the last section. 

We will now examine the influence of the number 
of samples, N, on the eStimation for open loop probability 
filter. The number of sampies has great influence on the 
estimation process because as the number of samples increase, 
the accuracy of the estimation increases. But, we have to 
look at the whole system, i.e.: 

- The dynamics of the process generating Z(t) 
when the dynamics of the process are rapid, 
we would like to take a minimun number of 
Samples. 

- A good estimation after N samples does not 
mean that we cannot get a better estimation 
with fewer than N samples (for example 
Nyaoy.s Lf we can filter the output of the 
probability filter by a Kalman-type filter 


we may get a better estimation. For example 
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N = 5 might be sometimes better than N = 1000 

if we place a Kalman-type filter after the 

PECbabiiatty filter. 
Thus, the information about the influence of the number of 
samples iS very important when we are going to design a 
whole system. The results of the test of the influence of 
the number of samples, are given in Table 5.2-6 and Figures 
5.2-9 and 5.2-10. From the results we see that the variances 
of the estimation errors decrease almost inversely (except 
for N = 2) with the number of samples. From the results it 
seems that the variances of the probability and likelihood 
filters decrease almost inversely with the number of samples. 
The likelihood filter (except at very small N) gives results 
for most of the cases better than the probability filter. 


d. Estimation of Parameter in Case of p.d.f. 
Not Known Exactly 


As was discussed in Appendix D, by having a refer- 
ence Signal we can estimate parameters evenif we have only 
an approximation for the p.d.f. (See equation D-12b.) 

In this section we would like to show how this 
concept works. We take at first a Gaussian process including 
a limiter, which iS very common in any practical measurement. 
We want to estimate the parameters with a minimum influence 
of the accuracy of the assumed p.d.f. We test this case with 
the probability filter with reference signals, and we compare 
it to the likelihood filter. The block diagram of the pro- 


Mess 1S given in Fig. 5.2-ll. 
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1.0 Jva r(ae2) 





0.1 
Be 
10 
ze 3 
10 
m---j=P (Probability filter) iN 
~ 4 @---j=1  (biklihood filter) . 
10 
N 
1 10 100 1000 10000 
Pig. 5-4-9. The variance of 9, as a function of the 


Sample's number. 
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1G 100 1000 10000 


muceese2- 10. ihe variance of om as a function of the 


sample's number. 


148 





ASFTWET YFTM SSed0ad sya jo weIbetp yooTq aug, - 


Tt-2°S °6tq 





Jaz Lue 





L=(Z)UWA 
0=(Z)3 








> UDEM 


SSad0ug ueissney 





x PU ULW 


nil (zy*4 


149 





For the probability filter the reference signal 
is taken to be a limited Gaussian signal and the assumed p.d.f. 
is GausSian. For likelihood filter the assumed p.d.f. is 
taken to be Gaussian also. 

We test first the case of a symmetrical limiter, 


| = . The mean of the output is not influ- 


me. , |Z 
max 


min 
enced by the limiter. The estimation of 35 1s influenced by 


Seen limiter. The bias error of 9, is strongly influenced 


2 
by the limiter. The following parameters are taken for this 


test. 


# of samples (N) = 100 ; 


feOoueensemble membexs = 100 ; 


From the first simulation (see Fig. 5.2-12) we can conclude 
that: 
Hines prOobapi lity filter 1s much better than 
the likelihood filter when ee SAE ewe): & 


Lm. me ciel t col tneoderi iter becomes a 
better filter than the probability filter. 
(11) The bias error of the likelihood filter 
depends very much on the knowledge of the 
beGet., witleyson Ehe probabliity filter 
SEEOmelt. 1s jalmost constant. 
The second simulation (with a non symmetrical limiter is 


performed), i.e., |Z aa . We choose Z| = -l, and 


aa max alg 


we change Laat The unsymmetrical limiter influences the 


Jitsy6, 
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variances and the means, so we have to check the bias error 
of estimating ay and 4: The results are given in Table 
7.2-/. The results show that the bias errors of the proba- 
bility filter are influenced weakly by the pv.d.f. while for 
the likelihood filter there is a great dependence upon the 
knowledge of the p.d.f. 

It should be mentioned that the likelihood filter used 
for comparison here is a "Straw man" in the sense that the 
meooa>oility filter is given the correct class of p.d.f's 
(indirectly through the reference signal), whereas the like- 
lihood filter is not. In this example, one could design a 
Meeeiinood filter using the correct class of p.d.f.'s that 
would probably outperform the probability filter. However, 
Smee pOolnt is that the probability filter is robust (in this 
example - further work is in order) with respect to errors 
in establishing the class of p.d.f.'s used for gradient 
computations when a reference Signal is available. In prac- 
tice, the choice between filter types would depend on whether 
the correct class of p.d.f.'s can be well approximated 
(likelihood filter) or not (probability filter) by simple 
analytic expressions. 

5. Conclusion of the Probability Filter Simulations 

After the simulations we conclude the following: 

a. It has been proven that the probability filter 

can be implemented quite easily. 

b. The performance of the probability filter is 


(in our example) not sensitive to errors in 
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estimating the form of the p.d.f. Ba) Oo when 
a reference signal is available. 

The performances of the filter when the p.d.f. 
is known exactly, is of the same order of magni- 
tude compared to the likelihood filter, but 
gives worse performance most of the time. 

In the target/chaff problem we do not know the 
exact p.d.f. because there are several types of 
chaff and targets. It 1s quite easy to provide 
a reference Signal for the probability filter 
WhGhn Wino thi ete the exact p.d:£. Thus the 
probability filter may be found to be the 
breaceucal Ssotution.» “further investigation is 


needed. 
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VI. SUMMARY AND CONCLUSION 


In this research we tried to solve a realistic and 
practical problem of separation of time signals when they 
have almost the same power spectrum. AS an example, we 
took the chaff-target problem in tracking radar where Doppler 
separation is impossible. We started with a discussion of 
the background of the problem and developed two possible 
solutions. 

a. The M.K.F. with antenna modification. This is one 
solution to the problem in the case of maximum 
information available (for example in the case of 
monopulse radar with phase detector). 

Beeeine probability filter. This is one possibility 
of the solution to the problem in case that minimum 
information is given (in the case of concical scan 
radar or lobe switching). 

To demonstrate our solution to the problem we have proven 
by simulations that the two main filters, the M.K.F. and the 
probability filters work as expected. 

The following items remain to be investigated as an 
extension of the work. 

meme ine complete solution in case of C.S. and L.S. 

b. The optimum number of samples (N) for the 
BEGoaoritty filter as a function of the dynamics 


of the signal process. 


ha) 





c. The probability filter extended to nonergodic 


processes. 
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APPENDIX A 


MoOoUtrP IED KALMAN FILTER FOR THE CASE THAT THE MEASUREMENT 
COMES FROM TWO DIFFERENT SOURCES 


i Lntroduction 

In many cases we have to filter data which is not 
the regular data assumed by the Kalman filter. The Kalman 
filter is well established and sometimes it is convenient 
to modify the problem such that it will fit the Kalman 
filter, so we will have a modified Kalman filter. In our 
case (See Chapter IV) such a problem arises. The problem 
can be described briefly as follows (refer to Fig. A-1l): 


we have a regular process described by (F-2) which gives 


us the output vector ¥ (K). ¥ (K) can be split into two 
vectors cy SL) (x) and 7 SU ae) ye 
v1) (x) 
¥(K) & |-------- 
vy §*) (x) 
where y (2) and y '4) are vectors of the same dimension. 


The vectors are transfered through a random process, 
which switches between them. After the switching block we 
get our measurement vector, Z(K). The vector Z (K) can be 


equal to one of the two possibilities: 


8) 


TK 
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bey) 
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Or: 


y ') (K) 
Z (K) ee naam 
frame's i 
The problem statement is: given the measurement Z(l), ..., 


Z(n) find the state vector X(n). 

Four papers deal with a Similar problem [26,27,28,29]. 
The first two papers cannot be implemented straightforward 
because both of them formulate the problem in a different 
way. The first one [26] deals with multi-target measure- 
ments giving information about the location of each target. 
The system must identify the type of the target and estimate 
its position which is a more general problem than ours, 
because we assume that we know the nature of the SESS 

The second and third references [27,28] deal with 
tracking in a cluttered environment. This could be adapted 
for our problem (by treating the clutter as chaff) but the 
assumption is that the undesirable returns occur completely 
at random which is not the case in our problem. There is 
no underlying dynamical process from which the returns are 
generated. Hence, no prediction from past data can be made 
On the location, nature or number of these returns at the 
next measurement time. 

The last reference [29] is close to our problem but 
Mses a different starting point. It does not give the 


optimum filter, which 1s known to be: 
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X(K|K) = ( sse6x a) dx 


— CO 


(for minimum variance estimates, see Appendix F). The 
vector 2 is all the measured data up to time K. The 
eea.t.'s e(x{2‘*)) are not Gaussian, which makes the evalua- 
ton Of X (K| K) very difficult to implement. It is extremely 
difficult to implement on a computer, since the p.d.f. is 
conditioned on all available measurement data. Thus [229] 
uses approximate p.d.f.'s. Our goal is to develop a filter 
Similar to the Kalman. In order to achieve this we use: 
- A selection block which transforms the measured data 
for input to a Kalman like filter. 
- The assumption that all the p.d.f.'s are Gaussian 
eeeke in) (29 ]) .- 
feacuseeel iG A-2- we try te reverse the switching 
in the data in the first block, by selecting one from the 
two possibilities of Z(K). Thus attermen:s block 1f our 
selection was correct we have the correct input to the 
Kalman filter. In a regular Kalman filter (R.K.F.) the 


input components gee 


(K) and Ze a) are not switched. 
Hence in our case if we made an error in the selection the 
output is biased. 

In our problem we have a positive probability to 


(1) (2) aS). he 


Select a wrong combination of 2 (K) and 2 
Baal modify the R.K.F. so that it will be able to handle 


ems error. The next block is a modified Kalman filter (M.K.F.). 
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Mis Lreads toua forum Like a Kalman filter but with 


modified gain and covariance matrices. In the R.K.F. we 


have two types of estimated states. 


One is X (K|K) which is called the update estimated 
State. This 1s the optimal estimation of X (K) taking 
into account all the observations so far (i.e., at times 
heece and including K). 

The second estimated vector is X(K|K=1) , which is the 
optimal estimation of the states given all the 
observations up to the last measurement (i.e., the 
observations at times up to and including K-1). 

This vector is called the predicted estimated states. 


The M.K.F. has similar types of estimated states, 


an update estimated state and a predicted estimated state. 
. ; 


Intuitively we can say that: 


If we trust our new measurement, we can use the 
update estimated state. 

If we do nottrust our new measurement, we have to 
use the predicted estimated state. 

If it is in between, the output must be weighted 
somehow between the two estimated states. 

The weighting depends on the probability of a 
SOrrect solution. 

The estimated states are not optimal, in the sense 
Of minimum variance of error, because of our 


assumption that the p.d.f.'s are Gaussian. 


i Gee 





2. Modified Kalman Filter 
Assumption 
(1) We have two separate, independent sources of 


States with a given dynamic: 


x) xt) = ox) gy eg ag cK) + WO (RK) 
(A-1) 
x‘*) (x41) = sy WSU ea (6 + A*?’ucR) + wt? 
(i1) The measurements are given by 
z(x) = u(x)[c'?? c'4)] }------ + V(K) 
e (A-la) 
IL Ok 
= L Ale.) orl 
Fo Jay sem 
L(K) = 
(A-lb) 
o 6m A, | 
- 0 = L. ate H., 
(111) The statisticof the random variables involved 


here are (assume all are Gaussian processes) : 


E(V(K)] = 0 ¥ K 

E(IV(K)V(J)] = R(K) d,, ¥% K 3 (A~2) 
wh) (x 

a) ==S=5 = 0 ¥ K 
BAP 


iGiZ 





E(W(K) *W" (K)] 
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A 
= Q(K) OL. AS 3 
= X 
~O 
COVAR Ss 
T 7 
E[X(0)-V'(K)] = 0 ¥ K 
ab 
ew Co) = 0) yk > 0 


i 


Egon (H-} = Prob (H,) = 5 


TL 
No correlation exists between H(k) at time K to H(j) at 
meme 7 for all K # j. j%In other words, there is no 

a priori information about the switching between 


Hy and Hi ° 
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For convenience we define: 


Fe) 1G) 
(i) gee 
Oo! 8) 
gst) ; 0 
bot | ----- F---- 
0 y 62) 
ye (A-2a) 
, & ae 
a 62) 
wi) (x) 
w(K) & | ------- 
w'?) (x) 
Solution 


We would like to have a sequential estimator 
because we like to estimate in real time with minimum delay. 
The non-sequential estimator (BayeSian approach) will be 
more accurate (i.e., with minimum variance error) but much 
more complicated, needs more storage, and has delayed output. 
For these reasons we will restrict ourselves to sequential 


estimators. Also, we will use a linear estimator. Because 
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the Kalman-like filter is linear, and we have to Store only 
the last information about our estimation, we will build 
our estimator in two steps - the first one will be a selec- 
tion block and the second one will be a modified Kalman 
filter. Because the Kalman filter is optimal for second 
order statistics and is a linear estimator, the only 
thing that we have to store is the previous values of the 
mmycOVarlance matrices and (i1) Predicted state. 

The modified Kalman filter has two main blocks 
(see Fig. A-2): 

a. The Selection Block. 

In this block we make the decision H = H, or 


H. - Let's define 


(a) Oy IG) 


~ ~~ 


i 


eS) 
(2) 


tk 
No 


By this definition we can write: 


See sla Hy is true 


iN 


(kK) = Ly 


---- Tf H is true 
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“TOFTTS UPUTeY PeTFTpow syi Jo wesbetp yootq ‘“7-¥ *bTa 


S19a1S psa ,.epsadd 


au}, jO XTJ1eW asueT4eADD 


a 





Sa1e1S pa}oapoud 


Cig 
UeWTeY We) 
pets TPoW) 
i vammmermesrtmeemanaes BuUTI4a} 1 ceeonnesc ana U0T}OeTAS a - 
¥ a EE OD nd 
V 


(1) 
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(see Fig. A-l). 


When the processing continues we have the predicted states 
(measured state) with the covariances of the error associated 
with the predicted states (measured state). Then by the 


likelihood ratio we can decide Hy On Ho - 


(Bae ee 
is = = (A-4) 
5 

£(z?) ota.) 


Qo ox 


Wc 
rh 


Y 


Thus the decision is: 


1 ah decide Hy 
(A-5) 


ie <a, decide H. 


By the knowledge of the covariance matrix, and assumed 
Gaussian distribution of the state, L can be computed. 
The predicted measurement is determined in the 


filter procedure: 


A(K) = the covariance error of the predicted 
measurement, Y(K|K-1) (A-6) 


G7) 





25 (K) 
2 ee. 2 page) Se ; the predicted measurement 
2) 
25 (K) 
By defining 
ri(K) §& ¥ (kK) - t, (x) 
~1 ~p ~ 
r. (Kk) = ¥, (X) = t. (kK) 
eS) 
Zia 
A ae 
cS i 
a‘) | 
Ae 
A oe 
es!) 
a! 
L.(K) = rT (k) A(R) +x, (K) 1 — ae ee, 


aus : 
(alr) 


(27) 2°) 1.) = (20) 4] a[-)/2expt- § 1, (K)} 


ia aimenasion Of Ehe measurement vector. 
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Substituting into Eq. (A-4) we get the decision law: 


L, (K) < L. (Kk) decide Hy 

VAS S)) 
L, (K) > L. (K) decide H. 
See the block diagram in Fig. A-3 for details. 

bee tne Bailtering Block 

After the decision about the measured state is 
made, we have now the problem of filtering the data. The 
Filtering could be straightforward if we have high confidence 
that our decision was correct, which is assumed in the 
regular Kalman filter. In our case we are not sure if our 
decision was correct, and so in the filtering process, we 
must take this into account. Because we cannot tell if our 
decision was correct or not we cannot reject or accept the 
data, we can only "weight" the data and depend on the 
computed probability that our data 1s correct or not. 

We have two possibilities: 

(i) A true data available, i.e., we decide 
the right decision, we will call this event D(K) - the 
data 1s desired, in sample K. 

(ii) A false data available, i.e., we decide 
the wrong decision, we will call this event F(K) - the data 
is false, in sample K. 


So we can write: 


JES, 





£(X(K)|Z(K)) = a(K)-£(X(K) D(K),2Z(K)] + [1-a(K)]£(X(K) |F(K) -2(K)) 
(A-9) 


where 


a (K) Prop decision K was desired |} 


Prob{D(K) |2(K) } (A-10) 


Assume that is Si Oe Then: 


a(K) = P{D(K)|Z(K)} = P{H, (K) |Z(K)} 


£{Z(K) |H, (K) }*P(H, (K)) 


Bene (Hy (Ky) Fe (z(K))H, (K) -PaH, (K) F 


(A= 1s) 
Observe that: 
SCH, (Qi =  wibs, <tr 
i D Z 
Then: 
(Ck) a er) (A-12) 
+ 6. (K) / 
where, in case of L. (K) > L, (K): 
=~. % £{Z(K) |H, (K) } _ eh) t? 
iL —  £{Z(K) 7H, (ky I 
and in case of L, (K) > L,(K), then: 
«wy 4 £{Z(K) |H, (K) } Gee 
ee SiZzCKa(K)) 


ab ye, 





Then we can conclude that: 


a(K) = spe (A~13) 
where: 
aie) & eee (A-13a) 
From Equation (A-9) we see that: 
X(K|K) ee Oe) ee (En X[K) | (K) , 2 (K) } 


(A-14) 


coo mac oe I) | Eu Zi) 


i.e, the new estimate of the state will be weighted with 

the results of a filter when assuming true data, and the 
predicted estimate of the state (assuming false data). The 
guestion now is, what kind of filter to use here? It seems 
reasonable to use a modified Kalman filter, because Kalman 
filter is the best linear sequential filter (when only the 
second order statistics are given). We use here the notation 
“modified Kalman filter" because it is not a regular Kalman 
filter in the sense that the covariance matrices and the 


gain are data dependent (because the output is a function 


eval 





of a). Furthermore, there is additional selection block 
which does not appear in R.K.F. 

The only way not to have extremely time con- 
Suming computations is to assume that the p.d.f. is 
Gaussian. The Kalman filter gives us the information that 
is used to evaluate an approximation of the p.d.f., assuming 
a Gaussian distribution. 

Because our gain and the covariance matrices 
are data dependent, we must compute them on line rather 
than off line. 

Let's summarize the results first and then 
derive them: 


a. The estimator is characterized by: 


A N AN 
Pata ee | Keb) iG Ae) [e4K) > CX(K)}K-1) | 
(=>) 
Bek had} = OX (K=1 | K-i) + AU(K-1) 
where ct is the input vector after the selection. 
b. The gain is given by: 
G,, (K) = a(K) G(R) (A-16) 
where 
G(K) = P(K|K-1) C° [CP (K|K-1) eaaet lo 


G(K) - is the regular Kalman filter gain. 





c. Update covariance error matrix, and predicted 


error covariance matrix: 


P(K+1|K) = 9%) (K[K) ot + Q (A-17) 


)(K|K) = 0” (K) P(K|K) + {1 - 0 (K) ]°P(K|K-1) 
+ 20(K) (1 - a(K)J{(I - G(K)C] o) (K-1|K-1) 67 


+ [I - G,(K)C1Q} 
where 


P (K| Kk) 


[I - G(K)C] P (K|K-1) 


P (K|K) the regular Kalman filter. 


lbw Pek 


_ a [Pr (K) Ly (K) |/2 


B (K) (A-18) 


a, (K), L» (K) are given by (A-7) 


e. Initial conditions: 


P(O|-1) = M (A-19) 





CO ath) Xo 
f. The predicted meaSured error covariance 
Matrix 1S given by: 


A(K) = CP(K[K-1) cl + R(X) (A-20) 


(See Fig. A-4 for details.) 


Explanation and Notation 


1. The differences between the Kalman filter 
and this filter are: 


-- The gains are data dependent, and so they 
must be computed online. 


-—- The error covariance matrix is not the 
same as in the Kalman filter. 


-- We do not use the measured vector Y but 
Ge Vector a2. i 


2. The Modified Gain - Gui 
In the regular Kalman filter the estimate 


equation of the state 1S given by: 


me (KR) = X(K|K-1) + G(X) [¢(K) - CX(K|K-1)] 
(A-21) 
X(K[K=1) = $X(K-1|K-1) + AU(K) 


where 


xR) cgi gy is the update estimate state in R.K.F. 





Po, 0 NOTTS 
IHL 





Butdwaqtty 


SIldLl 





caer adnpad0dg 





qUaWwUNsay 


OM} 3u} 





Ne 


out 
wWuo4 








| Se Eee 


T~ 


cP. 


V 


ruopaz0ud Butvaqytty ayy Jo aweuberp yootq ayy ow “Ord 


| 


1 a a ee, 


UNdAIOUd ONTYSL WS 





Ja} [td 





eo 


URWT ey} 





PeTLIPON | ‘Juorzoaqes 





— —————————— Ons 





In our case, when we compute the state X(K|K) we must take 
into account both the updated states assuming that our 
measurement is correct, and the predicted states assuming 


that our measured vector is not correct (A-14,15), i.e.: 


X(K[K) = a (kK) XS) (KEK) + [1 = a (K)]X(K|K-1) 
(A-22) 
ee eee Kec) Guk) [6 (K) =~ Cx(K|K-1)] 
Thus we see that the modified gain is Gy (®) = a(K)G(K). 
3. Covariance matrix of the error: 
The covariance matrix of the error is 
defined by: 
A Ak 
DOSS) Sia [e| 18) OG O| eg) (A=23) 
And: 
A 4 
red ee ea |) XC) 
= a(K)X(K|K) + [1 - a(K)]X(K|K-1) - X(K) 


Substituing e into }) we get: 
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(KK) = a7 (K)P(K[K) + [1 -a(K)]2P(K|K-1) 


+ 20 (K) [1 - a (K) JE(e(K[K)e™ (K|K-1)] 





where P(K|K), P(K|K-1) are the regular Kalman filter error 
covariance matrices. Now we have to evaluate the estimation 
of 
T 
e(K|K)e (K|K-1) ; 
It can be shown that: [35,38]: 
e(K|[K) = (I-G,(K)C]¢e(K-1|K-1) +G, (K)V(K) 


el: = G, (K)C] W(K-L) 


© (| oxoat) SOR IE het AS) 


maus : 


e(K|K)e> (K|K-1) { [I-G, (K) C] ge (K-1| K-1) +G, (K) V(K) -G(K) V(K) 


ie e yc) a (Kaaba de (K-1 | Ko Bye W(Kad) } 


= [I-G,,(K)C] $e (K-1|K-1)e° (K-1|K-1)9" 


Au sh 
+ G, (K)V(K) e> (K-1|K-1)® "= [I-G,,(K) C]W(K-1)e (K-1|K-1) 6 


[I-G,,(K) C ] $e (K=1 |K-1) W (K-1)-G, (K) ¥(K) W (K-1) 


+ [I - G (K)C]W(K-1)W? (K-1). 


—_— eee, —e 








Using the properties: 


E{e(K-1|K-1) v’(K)} = 0 
1 
ee eaiecK) {= 0 
a 
E{W(k-1) i) J 
ve get 
E{e(K|K)e (K|K-1)} = [I-G_(K)C] @] (K-1| K-1) 


+ [I - Gy(k)C]Q 
Substituting the result into the equation for }) we get: 


(K]K) = a7 (K)P(K|K) + [1 - o(K)]7P(K]K-1) 
+ 2a(K) (1 - a(K)]{[I-Gy(K)C]¢) (K-1|K~1) 6" 


#5 (LSE SIO = 
M 


4. The predicted measured error covariance 


Bisex (A (K)): 
OA * S it 
Coe sul (kK) Ko1) > Y(K)}1Y(K|/K=1) — ¥(K)]"} 


¥(K|K-1) = CX(K|K-1) (A-24) 


lee 2 





¥(K[K-1) - ¥(K) = C{X(K|K-1) - X(K) - V(x) 


From Eq. (A-24) we get 


a 


Chg eee OmneroK—) e at RK) 





These equations are Simulated and tested as discussed 


in Chapter IV. 





APPENDIX B 


THE MULTIVARIABLE NORMAL DISTRIBUTION 


In this dissertation the multivariable normal distribu- 
tion appears in many cases. As with the M.K.F., this 


appendix summarizes the equations of such a process. 


Let's denote the random variables Kproceer X, by the 
column vector: 
pt 
X = : (B-1) 
x, | 
The covariance matrix A is defined as: 
T11 Jia cecceee One 
Pt 
A = (B=-2) 
are orn 
al 
Where the elements O55 are given by: 
= 7 aCe B=-3 
T5 5 eee Ce eae oe) |) (B~3) 


JeSok 





The matrix 1S symmetric because 


ae = O54 (B-3) 


And the vector u (the vector mean of x) by: 


Te 


ze! . = (B-5) 
¥ on | 


It can be shown that, for the n-varliate case, the joint 


ee 


density £ (x) 1S given by: 
£(x) = (20) | a|7?/*exp{- S(x-u) "At (x-u)} (B68) 


where 


ee 


| A det {A} US 7) 


Also, we can write in matrix notation that: 


A = B{ (X-U) (X-U) *} (B-8) 


CG 
i 


E{X} 


Ikis) 2 





AS a matter of interest, we note that the n-variable 


normal moment generating function can be written as 


Ce) exp{U't + tAt/2} (B-9) 
where: 
Ss, 
t a (B-10) 
5 
n_ 
T 
wu) Efexp(t X) } (B-11) 


And by direct differentiation we can find that: 


atbt+... 
oe a eee ae Sai) 


pe eat 


b 
Z 


a 

E{X, x 

The results are for the four moments, in four dimensional 
case: 

E(X,X,X3X,] = E[X,X,]E[X,X,] + E[X,X,]E[X,X,] 

(B-13) 


ts E[X,X,]E[X,X,] 2E[X, JE[X,]E[X,EIX,] 


G;3 


One of the interesting and useful properties of normal 
random variables is that they are invariant to linear 
transformations. In other words, a linear transformation 


of a normal random variable is a normal random variation 


EOO. 
The result is that for the transformation 
Y = AX (B-14) 
where: 
al 
Y = : 
' Y ij 
n 
A =nxXm matrix 
is 
A s 
Uy =e al = A (B-15) 
A ay 
My —EGOUAR LY |.) = |ane (B-16) 


Let's take a special case but very useful for bivariate 


normal distribution with: 








VARIX,] = VARIX,] = 0° Gey 
Bie) om p 
leeeZ 
then: 
iL p 
A = a? 
QO 1 
Em@emwe Can write the p.d.f. by: 
‘ 2 y) 
z iL 5 Xy +X. +20X,X, 
fy (XX) - 5 18) a2) |= <a (B-18) 
270 AO 0.7) 
The p.d.f. of one variable given by 
2 
Ey. (X) = te epey | |S as] (B-19) 
ie 2 JACI 
ne 
fgemecOndLtional p.d.f. is: 
eee 2) i) - thx) (B-20) 


Thus 


185 





ae 24x 2 42x Xx 


aie eee aaa ie 


iL 
ie ae.) a 
Den |X Zero D 2 
aut \/ 200% (1-p") A ae 


] 


(OX, +X 


\? 
= etait a exp [- 1 2 


— OS] 
2 2 
fono* (1-97) 20 ({(l=o0 ) 


(S25) 





APPENDIX C 


TRANSFORMATION OF RANDOM VARIABLES 


Consider the problem that arises when a random variable 
x 1S transformed to a new random variable y through a 
functional relationship y = h(x). The usual problem is, 
megeeecene distribution of y, given the distribution of x. 


Bema monotonic function, it can be shown that [22], [1]: 


ep) SOY) S| (C-1) 
Seeby writing: 

x = h -(y) (G20) 
Then: 

£(¥) = fyfh “(¥))|$pth7*+ (x) 1 | (C-3) 


This only applies when h + SScrsecmencmeate ts Continuously 
differentiable function of x. For this case, we must 


eease fLrom the basic idea which is: 


Enesco Dye = (Set of x Corresponding to Y < b) (C-4) 


Tes 








Meeeexample, the set of x corresponding to Y < b is the 


intervals 


(—9,x,), (X5 7X) 


So by the definition of £Y (x) we can write: 


he 


Fi(b) = Jf£,(X)dx (C-5) 


where 


rT = fall the values of x for which Y < b}. 
A very useful example of transformation is the square 
law transformation: 


y = x ie oy) 


nen : 
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7b 





Os) ee CO ab) > 40 (C-7) 
-vb 
differentiating with respect to b, we have: 
- f(b ey i= aD) 
fy (b) = 
2Vb 
Replacing b by y we get: 
oe faeces VY) 
fy (Y) Oe a (E=3)) 
2vY¥ 
Now let's take the normal Gaussian density; then we get: 
f£y(¥) = —=£ [exp(-¥/2) + exp(-¥/2)] 
2V27TY 
Or 
f(y) = (2ny)7/* exp[-y/2] (C-9) 
Now, let's take the case of two variables: 
ieee We 5” 
(C-10) 
ee 5) 
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It can be shown that if Vi and Y5 are related to x] and 


Xo through continuously differentiable transforms then 





al 122): 
eo es «| 4) | (C-11) 
where: 
3X, 3X, 
oY 5 oY. 
J(x|y) = det (2) 
aX, 3X. 
ome oe 


J is called the Jacobian of the transformation. Where the 
transformations are not one to one, we must start from the 


definition: 
eB U55 gi)! j J (Ky Xo) dx, dx, (C-13) 


r = {set of x such that ue b 


side ra) 


and y, < b,} 


Now, by differentiating F., with respect to by and b. and 


Ne 
setting by aa b, = Yo, we get the proper p.d.f. 
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REE ENDE © 


SEPARATION OF SIGNALS BY PROBABILITY FILTER 


meee te COdUCt1iOn 

The general problem of separation of two signals is 
well discussed in the literature and the most familiar 
are Kalman and Weiner filters. We, in general, call one 
Signal "noise" and the second signal is called simply 
"Signal", and the problem is to estimate the signal when 
we have a noisy measurement. The most discussed and 
developed method of separation is for linear models, l.e., 


in control theory we have a general discrete model: 


eee Doty ee Viaje T Wid) 
(pe) 


2(i) = C X(i) + D V(A) 


When Xx is the state vector, Z is the measured vector, W 
and V are noises, U is known control inoue Deeae, 65°C, D 
are matrices (may be a function of time). 

A problem arises when the system is not linear. There 
are many possibilities to solve this problem, and the most 
Widely used are the extended Kalman filter, maximum likeli- 
hood filter, and the Bayesian approach. Here we use a new 
method which will be called the "probability filter." The 


advantage of this filter is that the knowledge required for 


ee a 





this filter is less than for the likelihood or Kalman 
filters, i.e., the only required knowledge is at most the 
a priori p.d.f. of the measured Signal (Sometimes even the 
knowledge of the p.d.f. might be non accurate). We will 
show that this method is an extension of the method known 
as "method of moments." In the next paragraph we will 
describe briefly the method of moments and after that we 


will introduce the new approach. 


2. Review of the Method of Moments 

The method of moments iS simple and intuitive. It is 
practical in the sense that it leads to a computationally 
reasonable estimator. This method can be established without 
E@emd priori knowledge of its p.d.f., but it Yequires a 
conditional p.d.f. on the obServation. In practice, it 
requires even less than this — only knowledge of the several 
first moments of the process. It yields an estimation that 
is not necessarily optimal in any sense. On the other hand, 
the estimate approaches the true value as the amount of 
data processed becomes infinite. So, on the one hand, it 
is sometimes slow in obtaining results because the amount 
of data required for the estimation is high, but, on the 
other hand, in many cases the result is guaranteed to estimate 
the true value. We conclude that, in general, this method 
Simply estimates the true value, but it is not an efficient 


estimator. 


oi 





Now let us introduce the basic concept of this method 
(according to Ref. [1] and [51]). The two references do 
not refer to the ambiguity problem which we have discussed 


above. Suppose that we have unknown vector parameter 


FS 
i 


[6,5 ca ey 8! (D-2) 


and we have sampled data from Z(t); 


I 


Z(t.) , j= 1,2, ..., m 
(D-3) 


Z 


a 


which depends on the vector parameter. Let's assume that 
the knowledge of the conditional moments from 1 to J of 


Z(t) is given, and exist, 1.e., 


1 ee 


Ot Biz, {6} = £,(6) , i=1,2, ..-, 4 (D-4) 


al 


Let's suppose now that we have J estimators of the cae 


moments, sO we can write: 


R 
lI 


l £, (9) 


(b=5)) 
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Now we have J equations with n unknown parameters, 85795, 

ay 8° In general J > n because the functions are not 
linear. 

If we have the same number of equations as unknowns 
it can lead to a non-unique solution. So for uniqueness 
we can use more equations than the number of the unknowns. 
Note: The problem of ambiguity can be also addressed by 
the Modified Kalman filter (see Appendix A). 
Bxample:  D-la (Reter to Fig. D-1.) 

> (8), with two 


we have two possibilities for the 


If we are given only two equations, £5 (8) and f 


parameters, 9, and 9 


Pil 
and g 2) 


il 
a) 


Weeror 6: 6 To decide what is the acceptable 


solution, we must have at least one more equation. With 
the addition of a third equation we accept one of the 
solutions. Note that the third equation does not lead to 
the same solution as the first two equations, because of 
the error in estimating the moments. BecauSe in many of 
the cases, the eStimation is an increaSing function of J, 
the last equation is used only for decisions and not for 
SStimation, i.e., our solution is: 

Ck), 


) = a, | choose 4 


es 


(2) 


(2) A 
Beit ,(8@ )-o choose 9 


The block diagram of the method of moments is given in 


Beg. D~-2. 
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In most of the cases the computational method of 
estimating the moments are simply by the definition of 


the moments, i.e., 


>) 
jt 
m1 9g 


eae (D-7) 


It is not difficult to generate examples in which the 
method of moments does not even lead to an acceptable 
estimator (see Ref. [{1]). In the following section we will 
see a new method which is a generalization of the method 


of moments. 


meee = Probability Filter" 


In the method of moments we saw that when one wants 
to estimate the parameters the first gee moments of the 
observed data is taken and we equate these to the J functions, 
the moments, and solve for the unknown parameters. Let's 
now suppose that J > ~. We choose the parameter in such 


a way that these infinite equations will lead to minimum 


error in some sense. For example: 
> [a = £.(8)]* ave mE (D-8) 
a ES) 1 
1=1 
when Ws & Oe In general, Ws 1s chosen as a decreasing 
function of i, for example w, = l/varta, - £,(6)}. The 


bo, 








method of moments is one example of D-8, i.e., by choosing: 


i 8 Pea 4 “nae! 

w. = €. > Q LP ie ae ee (D-9) 
a. 1 

w. = Q > a 
BL 


We get exactly the method of moments discussed before. Now, 
let's understand what we mean by Equation D-8. If we take 
W; > QO *1, then we are trying in some sense to minimize 
the difference between the conditional p.d.f. and the 


measured p.d.f. That is, we can modify Equation D-8 by: 


(2) || > min (D-10) 


where 


A 


when £, (4) denotes the estimation of the p.d.f. of the data. 
(We will see later that for the probability filter we do 
not even have to estimate the p.d.f.) In the last equation 


we try to choose @ in such a way that the norm of the differ- 


ence between the two p.d.f.'s 1S minimized. 


nos 








Een le D-2. (Refer to Fig. D=-3.) 
Let's assume that we have only one parameter to esti- 


mate; from the figure we choose 4 sigvel iplone 05 because 5 


gives a theoretical p.d.f. which is closer to the measured 
emer. than 4: 


Let us now choose the norm: 


fo ©) 


I= 5 f [£,(2|8) - £,(2)]° 42 > min (D-11) 


= © 


wn 


That is, we choose 9 in such a way that it minimizes the 
difference between the two p.d.f.'S in a mean Square sense. 
Now, let us take the derivative of I with respect to 

6, and set it equal to zero to find the minimum, i.e., a 


~ 


necessary condition for minimization is: 


= i - af, (Z| 8) 
0 = — = Je Ry CAI) a) eas 
36. ~ 36. 
SL =O 28 
" 9¢ (2/8) " ais (A eles 
ef 2 tale) az - ff En) 
96 96. ~ 30. 
aL -—cC al =O 1 
(ors 
a o£, (Z| 8) m o£, (2| 8) 
= g,(8) - —————— = g,(8) - —-—— (D-12a) 
30. i 30. 
als AE 


nog 








"€- PANHT A 


sie (laa 


200 





where the symbol "~~" means time average. The time average 


can be approximated by: 


a i 
df, (Z| 9) 1 2 9f,(2(t) |8) 
= = | er. Fe (D-12b) 
385 0 a6 5 


In our problem we try to estimate the amplitudes 
received at each receiver (which is porportional to the 
RCS of the target, see Chapter IV.A). Since we restrict 
ourselves to surface targets (which means low dynamic), 
the RCS of each target (ship and chaff) can be assumed to 
be stationary compared to the missile's flight time. 
Hence we can use "T", the time interval at the last 
equation (D-12b), small enough compared to the missile's 
merohit time. 


Eq. (D-12,D-12a) follow because: 


- ~ 3£, (2|8) ; 
g.(8) = f[ —=x— £,(2|6) az 
: 36 : 


Le i 


is a known function of 6 and can be computed a priori. 


Also 
2 | BK 
- £2) azZ 
piel 
— OO L 
Due AI) 
is the mean value of the function ———— with respect to 
00. 
aL 


Z0 1 


OO 





Z. Because we are aSSuming ergodic processes, the ensemble 


mean and the time mean are the same. So we can write: 
es a 
i o£, [2(t) | 9] 
g, (8) —- —————- = 0 (D-13) 
~ 3G. 


Because we have n parameters, we have n equations. 
Note that also in this case as in the method of moments, 
the equations are not linear and we may get a non-unique solu- 
tion. Thus we may have to add more equations to the set of 
equations D-13; for example, the equations of moments as 
discussed previously. Alternately we propose to use the 
M.K.F. (Appendix A) with all the possible solutions. 

A block diagram of the probability filter 1S given in 
frame D-4. The input Signal Z(t) 1s substituted into 

9£, (Z(t) |8) 


n~ 


08. 
a 


ama the function is averaged by the block following. Also, 
g (8) 1s produced. If the difference between the two is not 
zero, the value of 3 is changed until the difference between 
the two approaches zero. The initial condition for the 
integrator can be determined by the method of moments or 


as in the Kalman filter, the statistical mean of the 


parameters. 
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Example D-3. 
In this example we introduce a simple case of estimating 
a parameter by the probability filter. We choose the 


Saeemential distribution. The p.d.f. is given by: 


18 eA) aS Z > 0 (D-14) 


We can get the equation for the probability filter: 


i= (74 | 15)) _ _ 
a ee (D-15) 
ele 
~ 8£,(2| 8) ‘ 
ieee Dc edz = + = (D-16) 
0 
Thus g, (9) has a very Simple expression, 1l1.e., 
g, (9) = ob 2s (D-17) 
Meus, the equation for the probability filter is, from 
Eq. (p=—13), 
WG AE) SB ) e (O46 iL) | aa (D-18) 


i=1 


The last equation is the probability filter equation. It 
is a non-linear equation and we must solve it numerically 


mor 6. 
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Until now we have seen the direct implementation of 
the filter, i.e., determine the probability filter equation 
and solve for the parameters. Direct implementation of 
the filter 1s only one possible way of implementation. 
There are other possible implementations which give the 
filter its uniqueness. We introduce two possibilities of 
using the filter. In many cases we do not have the exact 
p.d.f. of the process but we have an approximation of it 
and a reference signal which represents the signal (e.g., 
by recording the signal previously). Wé give here two 
ways of using the probability filter when the exact p.d.f. 
1s not known. 

(1) The solution when only the approximate p.d.f. is 
given and a reference signal is available. 
This filter gives us the opportunity to estimate parameters 
even 1f our information about the data is not precisely 
given; that is, if we have a record of this kind of data, 
we can compare the two systems very easily. 


Let's return to Eq. D-12b. We have: 


= ae _(2|8) ee (218) - 
<< = (922225; (2/8) az of Bee 7B) ole 
Oe, - 38, an 38 

Sec 

= Q (D-19) 
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Let us assume that £ (2/8) is only an approximation of the 
Beoee., SO we Cannot “trust” this given p.d.f. completely. 
Like in the Kalman filter, we asSume an approximate model 
of the system, and by closing the loop we guarantee that 
the estimator behaves quite well even if our model does 
not precisely describe the real model. 

Equation D-19 gives us the possibility to close the 
loop aS in a Kalman filter. Eq. D-19 has two terms. The 
meeeemonec 1S computed “off” line if the p.d.f. £2 (2/8) is 
completely known. If the p.d.f. 1S unknown exactly, then 
the technique of comparing two terms is the best that we 
can do. We see that the two terms are almost the same. 


Assume we have a reference Signal XK with the p.d.f.: 


eae = alc) (D-19a) 
aus : 
A ae ae 
df, (2/8) ys df, (X|8) 
- £,(Z|8)dz = - 
96 ~ 36. 


Thus we can write: 


aes eas i Re 
ames) | AF (218) et pote) 
a8. a8. 00. 
1 al 1 
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The first term comes from the reference system, so by 


computing the two terms we can assume that when the 


w~ 


difference between the two terms becomes zero, § is close 


@emene true ¢. 


Figure D-5 demonstrates the idea. We measure the data, 


“~ 


Meee), transform the data to get the function h, (2,9), where 


nw) = — gui) (D=20) 
a9. 


‘ 
and g(w| 9) is an approximation of £(w|9) that has been 
chosen for its analytic form. 

In the original form of the probability filter equation 
(D-12a), we get an error when £,(¥/8) is not known exactly. 
Thus by implementing the probability filter equation (D-12a) 
we get a biased estimator. On the other hand if we have 

a reference signal, in case that f, is not known exactly, 
we Get: 


OO” ee a ee 


h, (X(r,9),9) - h, (2(t,9),9) = e; (Dae a) 
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For this case when 9 = 86: 
eq = Gln, (xX(t,9) 79) - h, (2(t,6),6)] (D-21b) 
= ferie (cqep mere cicmees 1) (2, o)t (28) dz 


This result follows because of the assumption that we have 
a reference Signal with the same p.d.f. as the signal, 
(pe2a) . 

Thus where the reference Signal has the same functional 
equation as the measured signal, we get an unbiased estimator, 
‘independent O38 h,(o,9), which is a function of the assumed 
p.d.£. However the feedback process can only be expected 
to converge when g (w| 8) is a reasonable approximation to 
£,(w| 8) 

In the regular "probability filter" we get a bias 
estimator when the p.d.f. is not known exactly. 

(ii) Indication of Change in a System 
Sometimes we like to check if parameters in a system have 
been changed. By taking a reference signal (from similar 
systems with nominal parameters or by recording the output 
from the system from the past), we can easily check if some 


parameters have been changed. 


ABS, 





Indication of change in the system is that: 


le.{| > ¢. (D=-23) 


where oy is a constant dependent upon the acceptable false 
alarm. The block diagram is given in Fig. D-6. 

By adding more calculations we might determine the 
value of the parameters that were changed, by simply 
eeeiying the probability filter. 

The probability filter implementation for the unresolved 
targets and the simulation of the probability filter are 


given in Chapter V. 


be be 











“A983 TTI AQTT Tqeqoid 


2943 Aq wajsdks ev ut ebueyo jo UOT}eOTpUT *9-q eanhTty 


ony 


: (2) 


WE}PSAS BOUGBLSL94 BY ST ZH WALSAS 
Wa1ShS paisa] eul ST | # welsks 





C oH 


WALSAS 


aBbueyud j0 


UO}yRIE Puy 





L ft 


7 | WALSAS 





A 


§ 
Vv 


(1) 9 





APPENDIX E 


THE MAXIMUM LIKELIHOOD ESTIMATOR 


A reasonable eStimate of a parameter is one which 

makes a given observation most likely, i.e., choose for 
the parameters 9, that value which maximizes the conditional 
p.d.f. £(2|@). The p.d.f. £(Z|@) is called the likelihood 
function. Because of the exponential character of many 
density functions, it is convenient to deal with the natural 
logarithm of the likelihood function - Inff(z[6)}. Clearly, 
because 217A) ()) ea ay lie et (A) 603 exists, so by choosing 9 to 
maximize In{f(zZ|6)} it maximized also £(Z|/9). A necessary 
condition that 6 maximize the likelihood function is: 

eae im (io) = 0 i = 1,2,.--.m (ea) 

a 

By setting all the m equations with m unknowns, one can 
solve for the unknown vector 9. It can be shown that in 
many cases the maximum likelihood function is a biased 
estimator, but the bias tends to zero [30]. Also the 
variance of the error tends to zero as numbers of samples 
increase. We can say that the likelihood estimator is 


asymptotically efficient. 
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APPENDIX F 
Powivotons POR SYSTEM DESCRIBED BY 
DIFFERENTIAL/DIFFERENCE EQUATION 
ifn TC RODUCTION 
In this dissertation we deal with estimators. The 
purpose of this appendix is to summarize the approaches 
given in the literature and establish the symbols and 
terminology used. There are two major approaches: the 
discrete and the continuous time. For our case, we deal 
with discrete time process. 
The difference equations for the discrete time and the 


measurement is given by: 


a = XK) RET KR) CCK) KR) WK) ) 
(a=) } 
bok) CU = © A (X(K),K) + V(K) 
where 
w,T,h - are known matrices 
W(K), V(K) - are independent sequences of indepen- 


dent random variables (not Gaussian 
in general) 


The only known case when the equations can be effectively 


solved (computed) exactly is the linear process plus 
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GausSian noise. This solution was obtained by Kalman [2] 
and Kalman and Bucy [3]. The known existence of this exact 
solution led to the approximate solution for the case where 
the plant and the measurements equation are almost linear 
and the noise processes are nearly GausSian [10]. But the 
Gaussian approximation and/or linearization are not adequate 
for many problems due to strong nonlinearity, as well as for 
non-GauSsian noise process. This has led to the development 
of many other kinds of approximations, for example, Gaussian 
Sum approximations [5], [6]. A good reference for many 
approximations can be found in the proceedings of the 6th 
nonlinear estimation theory symposium [7]. 

From equation (F-1) we see that the frequent problem 
is: 

The measurement noise and the plant noise are 

additive, with at most a state dependent multiplier, when 


@aeenmolse 1S in the form of Brownian motion. 


ue LINEAR SYSTEM — KALMAN ESTIMATOR 
Let's assume that our system described by a linear 


difference equation which is given by [3l]: 
Pa) eee 0K) AUK) + WK) 
(EZ) 


Y(K) = CX({K) + V{K) 


where 
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X(K) - are the states of the system at time K 


Di) ence Genero m anput to the system at time 


K (it iS a known, deterministic) 


W(K) =- the random forcing input 
V(K) - the measurement noise 
Y (UK) - the measured vector 


The statistical properties of the variables are given by: 


(1) The measurement noise: 


E{V(K)] = 0 K > 0 
eee 
E[V(K)v7(3)] = R6,- k,j > 0 
= oe ea ‘ — 
(2) The initial state: 
E(X(0)] = X, 
(b=) 
GOVAR WX(0)] = M 
eee ine forcing random process: 
15 ST aT 
(F-5) 
B(W(K)W'(j)] = Q6,, k,j > 0 


ES 





(4) The correlation between the variables: 
(a) The measurement noise and the initial states 


are uncorrelated: 


E(X(0)V'(K)] = 0 K 


|v 
7) 


(F-6) 


thie Gangom rosncing Input and the initial 


state are uncorrelated: 


E(X(O)W'(K)] = 0 kK 


lv 
7) 


(Eee) 


(c) The forcing random process and the measurement 


noise are uncorrelated: 
5 Lana ; 
EK) vtj) | a= 0 Ka 2 0 (B= Ss) 
The requirements for the filter are: 
(a) The estimator is characterized by linear 
relationship. 
(b>) Unbiased estimator. 
(c) The estimator minimizes the trace of the 
Sacer covariance Matrix, 1.e.: 


Trace {P(K K)} > min OPS) 


where 


Zr 





iieeresults of the Kalman estimator are given by [ 31] 


where: 


Notes: 


eo) 


- 


P (K|K) COVAR{X (K|K) = 0K) | 


X(K|K) is the estimated state depending on the 


input data up to and including sample K. 


K(K|K) = X(K[K-1) + G(K)[¥(K) - CX(K|K-1)] 
X(K[K-1) Ose ee (K=1) (F-10) 
X(0|-1) ee eo) 

G(K) = P (K|K-1)C? [cP (K|K-1)c?+R] 7+ 

P(K|K) = []- G(K)C] P(K|K-1) 

P(K+1|K) = oP(K/K)A’ + Q 

P(O|-1) = M 


The Kalman filter is the optimum linear estimator 
for the given knowledge of the second order 


statistics only. 


Pee] 





(11) There are other estimators, like nonlinear, which 
are Superior to the Kalman filter for the general 
case. 

(i1i1) If all the variables are Gaussian, the optimal 
minimum variance estimator turns out to be the 


Kalman filter. 


3. THE BAYESIAN APPROACH 
The BayeSian approach simply takes the expected value 
of the state, given all the measured data, l.e.: 


X(K|K) = E{X(K)|¥(0), ..-, ¥(K)} (eka t) 


Tf we define: 


RK 
A 
es 


i (0) eee LK) | 


Then (F-11) can be formed as: 


X(K[K) = ELx(K) [¥ 
If we have the p.d.f. of the state given re we can 
find the expected value of X, i.e.: 
X(KIK) = f x(x) £0x(x) |y¥"%)] axcK) GoD) 


LX] 
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Semcne problem is to find the p.d.f. 


ies £1x (x) [y 872) ye Ly (x) [x (x) 
etxcm [y= ———=—_ 
fe (SSCS) | se | 
etx(x) |x 2)) = ff eex(n-2) fy“8 1 y etx cK) [x(K-1) 14x (K-1) 
[x] 
ety(xy[y"")) =f etxcey fy ye ry cK) [x (x) ax) 
(X] 
etx(a)|¥°")) & ffx (0) 


Those equations, theoretically, give the complete solution 
of the estimation. The problem is that the integration 
Samenot be carried out in closed form. Thus, in most of 


the practical cases this solution is not used. 
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